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ABSTRACT 

A small fraction of core collapse supernovae (SNe) show evidence that the outgoing 
blast wave has encountered a substantial mass ~ 1 — IOA/q of circumstellar matter 
(CSM) at radii ~ 10^ — 10^ AU, much more than can nominally be explained by 
pre-explosion stellar winds. In extreme cases this interaction may power the most 
luminous, optically-energetic SNe yet discovered. Interpretations for the origin of the 
CSM have thus far centered on explosive eruptions from the star just ~ years— decades 
prior to the core collapse. Here we consider an alternative possibility that the inferred 
CSM is a relic disk left over from stellar birth. We investigate this hypothesis by 
calculating the evolution of proto-stellar disks around massive stars following their 
early embedded phase using a self-similar accretion model. We identify a brief initial 
gravitationally-unstable ( "gravito-turbulent" ) phase, followed by a much longer period 
of irradiation-supported accretion during which less effective non-gravitational forms 
of angular momentum transport dominate. Although external influences, such as the 
presence of a wide binary companion, may preclude disk survival in many systems, 
we find that massive 1 — IOMq) disks can preferentially survive around the most 
massive stars. Reasons for this perhaps counter-intuitive result include (1) the shorter 
stellar lifetimes and (2) large photo-evaporation radii (~ 10"^ AU) of very massive 
stars; (3) suppression of the magneto-rotational instability due to the shielding from 
external sources of ionization; and (4) relative invulnerability of massive disks to lower 
mass stellar collisions and luminous blue variable eruptions. Furthermore, disks with 
radii ^ 10^ — 10'^ AU are picked out by the physics of the embedded stage of accretion 
and the requisite conditions for subsequent disk survival. The disk mass, radius, and 
scale- height at core collapse typically result in a 10 per cent efficiency for converting 
the kinetic energy of the exploding star into radiation, potentially producing a total 
electromagnetic output ~ 10^° — 10^^ ergs. We identify two regimes of disk-supernova 
interaction, which are distinguished by whether the shocked disk radiates its thermal 
energy before being engulfed by the expanding SN ejecta. This dichotomy may explain 
the difference between very luminous supernova which show narrow H line emission 
and those which show no direct evidence for hydrogen-rich CSM interaction. Because 
very luminous SNe are rare, testing the relic disk model requires constraining the 
presence of long-lived disks around a small fraction of very massive stars. 

Key words: stars: winds, outflows, formation; circumstellar matter; supernova: gen- 
eral; accretion: accretion disks; infrared: stars 



1 INTRODUCTION 

Sensitive, wide-field transient surveys are revolutionizing 
our understanding of the breadth of time-dependent as- 
trophysical phenomena. The discovery of extremely sub- 
luminous, rapidly-evolving Type I supernovae (SNe) (e.g. 
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Foley et all l2009l: IValenti et all l2009l : iPerets et al] l2009l : 



Poznanski et al.l 20101 ) may, for instance, indicate a wider 
diversity of white dwarf-related explosions than 



ously anticipated (|Bildsten et al.ll2007f : iMetzger et al.ll2009l : 
iShen et al., ,2oTol '). In the opposite extreme, a number of 
core-collapse SNe have recently been discovered that are 
unusually bright and/or optically-energetic. These events, 
collectively known as very luminous SNe (VLSNe), include 
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SN 2006Ky jOfek et all l2007l: ISmith et al l l2007f ). 2006tf 
("Smith ct al. 2008'). 2005ap ("Quimb v et all boO ?'). 2003ma 
(Rest ct al._ 2009), 2007bi ( Gal- Yam et alJlToOQ ) . 2008fz 
dDrake et alj|200<l. 2008es (iGezari et al.ll2009l : iMiller et al] 
I2OO9I ) and 2 008iv jMiller et al."2010'). The mysterious event 
SCP 06F6 (jBarbarv et all 12009) and the luminous, high- 
redsh ift transients discove red by the Palomar Transient Fac- 
tory (jOuimbv et al.ll2009l ) may be related phenomena. VL- 
SNe radiate a total energy > 10^" — 10^^ ergs and can 
reach peak absolute magnitudes approaching My ~ —23 
(e.g. 2005ap). Though they represent only a small fraction 
of massive stellar deaths, VLSNe are thus particularly con- 
spicuous because they can remain bright for months— years 
and are detectable even in the distant universe (jCooke et al.l 
l2009l : lQuimbv et al.|[2009l ). 

The optical emission from most Type I and some Type 
II SNe is powered by the radioactive decay sequence ^^Ni 
— > ^^Co — > ^^Fe, suggesting that VLSNe may result from an 
anomalously high Ni yield (e.g. at le ast several solar masse s 
are required to explain SN 2006gy; ^Ka wabata et al.ll2009l ). 
Though such a large ^^Ni mass is difficult to produce in 
a normal core collapse explosion, it is a natural byprod- 
uct of theoretically predicted "pair-ins tability" SNe (PI- 
SNe: lBarkat et al.l 19671 : iBond et al.lll984l ). Although PI-SNe 
are a plausib le origin for some event s including SN 2006g y 
l|Smith et al.l l2007) and SN 2007bi (|Gal-Yam et all 120091 ) . 
they are an unlikely explanation for all VLSNe since PI-SNe 
are expected to occur only from extremely m assive stars 
with low metallicities (|Heger fc WooslevI |2002| ). as would 
form primarily in the early universe. 

Another way to power an anomalously bright SN light 
curve is to increase the explosion's prompt radiative effi- 
ciency; indeed, th e optical output of most core collapse SNe 
(~ lO**^ ergs; e.g. iBersten fc Hamuv , 2009) represents only 
a small fraction of the total > 10^^ ergs of kinetic energy 
that is generally available. Low radiative efficiency is typi- 
cal because the stellar progenitors of core-collapse SNe are 
relatively compact: radii for Wolf-Rayet and red giants are 
R ^ 0.01 and <: 10 AU, respectively. Much of the thermal 
energy generated by the shocked stellar envelope is therefore 
lost to PdV work before the ejecta expands sufficiently to 
become transparent. 

The radiative efficiency can be enhanced, however, 
if the kinetic energy of the ejecta is thermalized at 
much larger radii (> 100 A U), in a region ab o ve or 
closer to the photosphere (e.g. [Falk fc Arnetl Il973l . 1 19771 : 
IChevalier fc FranssonI 1 1994 [Smi th & M cCravl l20o'7f ). This 
can occur via shock interaction with a dense circumstel- 
lar medium (GSM). Indeed, th e Type Iln class of SNe (e.g. 
ISchlegel|[l990l : lFilippenkolll997l ). which are characterized in 
part by narrow H emission lines suggestive of slower mate- 
rial, are generally believed to be powered by GSM interac- 
tion (e.g. lGhugai et al.ll2004 ). Although many VLSNe are in- 
deed classified as Type Iln, others show little or no evidence 
for G SM interaction (such as SN 2005ap; iQuimbv et al.l 
|2007^ or hydrogen in their spectra. Such behavior may nev- 
ertheless be consistent with a GSM-powered luminosity if 
the GSM is especially massive and opaque, since this would 
delay t he time required for th e shocked thermal energy to 
escape (jSmith fc McGravl[2007l ) and for the SN to begin dis - 
playing more typical Iln properties (e.g. ISmith et alj f l2008l ). 
We describe a specific example of this effect in 8^5.21 



In order to e xplain luminous Iln SNe ranging from 
events like 2008iy d Miller et allbOlOl ) to extreme cases such 
as SN 2006gy via GSM interaction requires the presence of 
~ 0.1 - 10A/q on radial scales ~ 100 - 1000 AU. If one 
assumes that the narrow H line widths (~ 10^ - 10^ km 
s~^) commonly observed in SN Iln represent the unshocked 
outflow speed of the same GSM responsible for powering 
the SN luminosity, the inferred mass-loss rates range from 



much high er than can be 
Gavlev et al.l 



~ 10~^A/q yr~^ up to ~ Mq yr 

suppli ed by standard, line - driven winds (e, , 

1995; lOwocki et al] |2004| : ISmith fc Owockil l2006t ). When 
combined with evidence for minima or "gaps" in the mass 
distr ibution between the star and the GSM in some cases 
(e.g. iGhugai fc Danzigeij|l994 iMiller et al.ll2010l ). these in- 
ferences have motivated the interpretation that the GSM 
is composed of shells of mat erial ejected .just years— decades 
before the SN explosion (e . g . iGhugai et al.ll20o3 : ISmith et al.l 
I2OO7I : iDessart et al]|2009l ) . 

Proposed explanations for these pre-SN outbursts 
thus far include luminous blue variable (LBV)-like erup- 
tions analogous to_the 19th century eruption of r;-Garinae 
(e.g. ISmith et al.ll2007l ) and pulsational pair i nstabilities (i.e. 
sub-e nergetic, non-terminal PI-SN analogs; IWooslev et al.l 
|2003). However, neither of these explanations is altogether 
satisfactory. If LBV eruptions do occur just prior to stel- 
lar death, this suggests that very massive stars can die in 
their LBV phase , nominally though t to occur prior to core 
He burning fe.g. iLanger et al. 1994), an d with their hydro- 
gen envelopes intact fe.g. ISmithll2008bl ). Pair- instabilities, 
on the other hand, occur only in stars with masses > 95Mq 
that retain their He envelope following H burning; this again 
likely requires low metallicity, an unusual circumstance in 
the present-day universe. 

It has been argued that both the hydrogen-rich deaths 
of very massive stars and modern-day PI-SNe become more 
likely if empirically-calibrated stellar mass-loss rate s have 
been system atically overestimated (|Gal-Yam et al.l l2007l : 
ISmith|[2008 d) due to the ins ensitivity of n'^ emissi on diagnos- 
tics to wind clumping (e.g. iFuUerton et al.ll200^ ). This does 
not, however, mitigate the fact that both LBV and pair in- 
stability models also appear to require fine-tuning: the delay 
between shell ejection and core collapse must be well-timed 
to produce a collision at the optimal radii (~ lO'^ — 10^ AU) 
for producing bright emission. Although this "coincidence" 
could in part result from a selection effect, the fact remains 
that the pre-SN eruption and core collapse must be prefer- 
entially synchronized to within a timescale ~ years— decades 
(which is orders of magnitude shorter than stellar trans- 
port or evolutionary timescales) or VLSNe would he even 
more rare than is observed. Although such a correlation 
could in principle arise from the time interval between pair- 
instabilities, this depends sensitively on uncertain details 
such as ho w long the core takes to contract following the 
first pulse (|Heger fc WooslevI I2OO2I : I Wooslev et al.ll2007l ). 

In this paper we consider an alternative origin for 
GSM interaction in luminous core collapse SNe: a gaseous 
proto-stellar disk left over from stellar birth. The life- 
times of the mo s t massive stars are only ~ 3 Myr (e.g. 
iBond et aLlll984l : [Maeder fc Mevnetlll987[ ). comparable to 
the observed lifet imes of disks around low-mass proto-stars 
(e.g. IStromlll995l ). Observations of moderately massive (~ 
IOMq) proto-stars typically suggest much shorter disk life- 
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time s fe.g. iNatta et al.]|2000l : iFuente et al.ll2006l : see, how- 
ever, Im 

anoi et al indicating that the "disk dispersal" 

processes at work following star f ormation may be more ef- 
fective for higher mass stars (e.g. iHollenbach et al.|[200ol ). 

VLSNe are, however, very rare and likely originate from 
a small subset of very massive stars. The PTF events discov- 
ered bv lQuimbv et all (^2009) represent only a tiny fraction 
- 10"'' of core collapse SNe. Although Type Iln SNe as a 
whole account for a few to ~ 10 per cent of core collapse 
events (|Cappellaro et al] 1 19991) . they are a highly inhomo- 
geneous population and the majority may represent a differ- 
ent phenomena than the most luminous events. If a fraction 
/ ~ 10~* — 10~^ of massive stars must retain a massive disk 
in order to be consistent with the rate of VLSNe, then only a 
few to several hundred relic disk systems should be present in 
Milky Way-type galaxies at any time. Testing this prediction 
is nontrivial because although a relic disk would likely be 
conspicuous, massive stars are necessarily very young and, 
hence, typically far away and often obscured; observations 
and their interpretation are thus particularly challenging 
due to limited sensitivity, angular, and spectral resolution 
(e.g. ICesaroni et al.l [20071 ). Although there is no evidence 
for long-lived massive disks around unobscured ZAMS O 
stars, relic disk systems may not appear as normal or unob- 
scured, and the local census of massiv e obscured s tars and 
their environments is incomplete (e.g. iWachter et al.ii201G) . 
Indeed, the rate of VLSNe is sufficiently low that if only a 
very small fraction of high mass stars retain a disk, it would 
have important observable implications. 

The present theoretical study can broadly be divided 
into two parts: a detailed analysis of (1) the conditions un- 
der which a massive proto-stellar disk can survive until the 
core collapse of its host star; (2) the interaction of the relic 
disk (with properties determined self-consistently) with the 
outgoing supernova ejecta and the resulting emission. In de- 
tail, the following sections are organized as follows. In 321we 
discuss the stages in the formation and evolution of massive 
proto-stellar accretion disks. Relying heavily on the results 
from in 331 we construct time-dependent proto-stellar 
disk models which focus on the isolated stages of evolution 
following the embedded accretion phase, neglecting for the 
moment all sources of disk mass loss except accretion. In 
33] we examine the susceptibility of the disk to destructive 
processes such as outflows and stellar collisions in order to 
determine whether and under what conditions disk survival 
is most probable. In 3S]we explore the interaction of the disk 
with the supernova ejecta following core collapse and its im- 
plications for the origin of luminous CSM-powered SNe. In 
we discuss the implications of our results for long-lived 
disks around massive stars ( jjG.lf) . Type Iln (and non-IIn) 
VLSNe (El) and "hybrid" Type I/IIn SNe (iJOJ. We con- 
clude in 33 



2 STAGES OF DISK EVOLUTION 

This section provides a discussion of the formation and 
evolution of proto-stellar disks around very massive stars. 
We begin in ^2.11 by describing the embedded phase, dur- 
ing which the proto-star grows by active accretion from its 
parent molecular core. During this stage, most of the disk 
is gravitationally unstable; the outer disk is susceptible to 



fragmention and the inner disk rapidly feeds gas (and, pos- 
sibly, bound companions) onto the central proto-star. Once 
infall shuts off, the embedded phase ends and the remain- 
ing disk viscously evolves in relative isolation (neglecting 
for the moment external destructive processes; 331 • This iso- 
lated evolution generally begins, as during the disk's embed- 
ded evolution, with a gravitationally-unstable phase, during 
which "gravito-turbulence" supports the disk against frag- 
mentation and supplies the angular momentum transport 
f ^2.2|l . This phase is short-lived, however, because as accre- 
tion depletes the mass of the disk, stellar irradiation becomes 
an increasingly important source of midplane pressure. The 
disk thus rapidly transitions into a gravitationally stable, 
irradiation-supported accretion phase ( i]2.3|) . From this point 
on, other generally less efficient processes, such as MHD tur- 
bulence under partially-ionized conditions, supply the angu- 
lar momentum transport. This typically results in a much 
slower disk evolution until stellar core collapse. 



2.1 Disk Formation and the Embedded Phase 

The formation of a centrifugally-supported disk may be a 
nearly ubiquitous feature of stellar birth due to the sub- 
stantial angular momentum of the parent molecular gas. 
For low mass proto-stars such as T-Tauri sta rs, evidence for 
disks is strong and sometimes explicit (e.g. iBertoutI Il989l : 
iBurrows et~aiTl 19961 : TSimon et al1l2000l ). Although observa- 
tions are generally more challengin g for high mass stars (e.g. 
iBeuther et al.|[2006l : ICesaroni et"al.. 2007 ). advances in sub- 
millimetre observations of massive protostars (^ IOMq) 
have revealed the presence of flattened structures an d, pos- 
sibly, Keplerian disks ( e.g. IChini et al.| [2004: Ccsaron i et al.l 
I2OO5I : iPatel" et ahlliooH ). In fact, disk accretion is probably 
crucial to the very process by which very massive stars form: 
a disk shields infal ling material and re-directs the stellar 
radiation field (e.g . ^Stabler e t al.|[200ol : lYorke fc Sonnhalterl 
l2002i : [Krumholz et a l. 2005), allowing accretion despite the 
stifling effects of radiation pressure o n gas and dust (e.g. 
lKahnlll974l : IWolflre fc Cassinellilll98"7l ). 

During the earliest, "embedded" phase in massive star 
formation, the protostar accretes most of its final stellar 
mass M* 10 — IOOMq from its progenitor molecular 
core on a timescale tacc 1 — 2 x 10^ years, correspond- 



ing to an accretion rate M , 

-1 



10" 



lO-'MfT^ 





yr - (e.g. iMcKee fc Taiil l2003l : iBaneriee fc Pudritzl |2007| ). 

Such a large accretion rate cannot in general be accommo- 
dated by local visco us torques due to e.g. hydrodynamical 
or MHD turbulence (|Balbus fc Hawlevlll998l ). However, once 
the mass of the disk grows to a substantial fraction of the 
mass of the central proto-star, the dis k becomes suscepti- 
ble to gravitational i nstabilities (GI; e.g. lTomlev et al"]|l99ll : 



Johnson fc Gammieli2003. ). These generally set in once the 



Toomrd (|l964l ) parameter 



Q 



A/* H 



(1) 



ttGE tvR^T. R 

decreases below a critical value Qo ~ 1- Here Cg is the adi- 
abatic sound speed, f2 — [GMt, / B?)^^^ is the orbital fre- 
quency (assuming a Keplerian potential), E is the surface 
density, K'h is the stellar mass, tvR^T, is a measure of the 
local disk mass at radius R, and H — Cs/Q is the disk scale- 
height, assuming vertical hydrostatic balance. 
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Numerical simulations show that the nonlinear devel- 
opment of G I depe nds on the thermodynamic properties of 
the disk (e.g. lGamm ic 2001; Lodato fc Rice 2004). If the disk 
cooling time tcooi ~ SCg/F exceeds the orbital period 
where F is the disk's outward vertical energy flux, the disk 
settles into a quasi-steady state of "gravito-turbulence." Dis- 
sipation of this Gl-induced turbulence heats the disk (thus 
raising Cs) until Q ^ Qo, i.e. marginal gravita t ional stabil- 
ity obtains (e.g. [Gammi"3l200ll: lRafikovil2005l . l2007l . l2009l : 
iMatzner fc Le vin"2005 l: iBolev et al.ll2006l) . The effective di- 
mensionless [Shakura fc SunvaevI ( 19731 ) "alpha" parameter 
correspondin g to the requisi te level of turbulent viscosity 
was found by Gammid (|200ll ') to be 



9(7 - l)mcooi 



(2) 



where 7 is the adiabatic index. 

When, on the other hand, tcooi < , gravitational in- 
stabilities cause the disk to fragment into bound substruc- 
tures. One way to interpret this result from equation ([2]) 
is th at the disk cann ot provide stress at the level a > 1 
(e.g. iRice et~ai]|2005l ). The critical conditions for fragmen- 
tation, viz. ftgi ^ 1 for Q = Qo ~ 1, can be trans- 
lated into a maximum accretion rate A/max ~ Sttj/E ~ 
Cg/G that the disk can accommodate without fragmenting, 
where v = ggiCsJ/ is t he effective kinematic viscosity (e.g. 
iMatzner fc Levinll2005l) 



Kratter fc Matzneij (|2006l '). hereafter KM06, examine 



the conditions under which disks around massive proto-stars 
fragment during the embedded phase. In particular they de- 
termine at what radii the disk can support accretion at the 
core-supplied in-fall rate (viz. Ale <, A/max), taking into ac- 
count the stabilizing effects of viscous turbulent heating and 
irradiation from the central star. KM06 conclude that for 
typical core accretion rates Mc ~ 10""^ — 10~'^Mq yr~^ the 
disk fragments outside of a critical radius -Rfrag, which ro- 
bustly lies in the range -Rfrag ~ 100 — 200 AU for the entire 
range in stellar mass M-^ ^ 10 — IOOMq. 

KM06 then compare Tifrag to the typical circularization 
radius -Rcirc , which is determined by the angular momentum 
of the infalling envelope. They estimate that -Rcirc increases 
from ~ 100 AU to ~ 500 AU as A/* increases from ~ 10 
to ~ lOOAfQ, but with significant scatter at all values of 
Af* (in contrast to -Rfrag, which is approximately constant 
with -M*). From this they conclude that the outer regions 
of proto-stellar disks around the majority of massive stars 
are susceptible to fragmentation, presumably resulting in 
the formation of one or more proto-stellar companions at 
R > Rks-s- 

The presence of a massive companion at ~ 10^ — 10'' 
AU would likely preclude the long-term survival of a proto- 
stellar disk. A fraction of massive cores will, however, also 
form stable disks with 7?circ < i?frag from the low end tail 
of the core angular momentum distributiouQ which is broad 
because the core is supported by stochastic, turbulent pro- 
cesses. Furthermore, even when fragmentation occurs, it is 
unlikely to completely suppress accretion. Fragmenting re- 



^ In fact, by not sharing their mass with binary companions, 
such stars may preferentially grow to become the most massive 
(KM06). 



gions may be sus ceptible to additional, global gravitationa l 
instabilities (e.g. I Adams et all 1 19891 : ILaughlin et ahlirog^ ). 
which would cause matter to rapidly accrete t o radii < -Rfrag 
(i.e. with an effective viscosity a ~ 1; e.g. iKratter et al.l 
I2OO8I ). Thus, even if a proto-stellar companion forms out- 
side -Rfrag, it may migrate inside -Rfrag during the embed- 
ded phase, resulting in a tight bina ry or even coalescence 
with the central s tar (e.g. KM06; iKrumholz et all l2007l : 
iKratter et al Indeed, massive stars are known to have 

high binary fractions (e.g. IPreibisch et al.|[200ll : lLadall2006l '). 
with apparent mass-ratio and semi-major axis distributions 
which suggest that the binary orbits were formed during 
the embedd ed phase with small separations ^ -Rfrag 



100 AU (e.g.lPinsonneault fc StanekH2006l : lApai et al.|[2007l : 



IKrumholz fc Thompsonll200it ) 



In the rest of this paper we focus on the subsequent 
post-embedded evolution of the disk, assuming that its mass 
is initially concentrated near a radius ~ 100 AU < -Rfrag. 
This choice is motivated by the competition between molec- 
ular core angular momentum (which favors large radii) 
and the requirement not to fragment (which requires radii 
< iifrag). We furthermore neglect the possible effects of stel- 
lar companions at large radii. In doing so, we are implic- 
itly assuming that either (1) the disk forms from relatively 
low angular momentum material such that no massive stel- 
lar companions form; or (2) companions do form, but they 
subsequently migrate inwards during the embedded phase, 
resulting in a single star or tight binary with negligible sub- 
sequent effect on the disk evolution at larger radii. 



2.2 Isolated Gravito- Turbulent Phase 

As argued in the previous section, a plausible initial condi- 
tion following the embedded stage of accretion is a relatively 
isolated disk with initial radius -R^.o ^ 100 AU < -Rfrag and 
mass Mdfi (to be determined below). Although by construc- 
tion the disk is stable to fragmentation at its initial radius 
(i.e. fitcooi > 1), it still resides in a marginally stable state 
of gravito-turbulence with Q ~ Qo ~ 1, which we now char- 
acterize. 

Using equation (ITJ, the midplane temperature T of a 
disk with Q = Qo at the radius Rd where the local disk 
mass (X E-R^ peaks is given by 



Fgi ~ ^ = A50KQIA2^M^ 



woiMa/O.lM.fRa 



(3) 



where Af* = lOOAf*,iooA/0, Rd = 100i?d,ioo AU, 
Att-R^E is the total disk mass, and A = 2A2 is a constant 
that relates total disk mass to the local properties at Rd and 
which takes the value A = Agi « 2 for a gravito-turbulent 
disk in steady state (see Appendix A and the discussion in 
>13.1|l . The sound speed is Cs = {jkT/ ^^^^ , where we have 
taken 7 — 7/5 and ^ ~ 2.34mH for molecular gas, and mn is 
the mass of a hydrogen atom. Likewise, the disk scaleheight 
at -R ~ Rd is given by 

^ ~0.05Qo(A/d/0.1Af.)^"' (4) 
R Rd 

The massive proto-stellar disks of interest are dusty and 
optically-thick, with a midplane optical depth r = kE/2, 
where k is the dust opacity. The radiative flux through the 
disk surface is F ^ ScrT^/Sr (assuming a constant heating 
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rate per unit mass) , and the cooling time is thus given by 
l|Kratter et al.ll2010h 



32(7 - l)aT4 
1180yrQo''AUioM-loo{Md/0.1M^y 



(5) 



where we have scaled k = IOkio cm g~ to a constant 
value that is representative of the limited opacity range 
(~ 3 — 16 cm g~ ) appropriate throughout t he temperature 
range 10^ ¥^ <^T <. K of present interest l|Semenov et alj 
|2003| ). The gravito-turbulent viscosity (eq. [5]) correspond- 
ing to the cooling rate in equation ((5]) is given by 



Qgi ^ 1.2 X 10-2QgAj''K-^A'4\(o(,(Md/0.1M* 
corresponding to a viscous timescale 



,5/2 



„(6) 



Q\Rd — Qo 



0.03Qo''^K~'''^M 



-1/6 pll/6 /, /, xl/4 /, rx 

*,100 -"-0,100 l'-/ '■vise, Oj ) l-LOJ 

(16) 



where we have used equation ((§]) for the initial disk mass. 

From the above, note that the disk accretes and spreads 
slowly in time {Md oc t~^^'^;Rd oc t^^'^) due to the rapid 
increase in the viscous time with decreasing disk mass 
(tvisc (X. M^^; eq. [7]). Also note that (fitcooi)"^ oc Ogi 
increases with time as t^'''^, which suggests that the disk 
may eventually become unstable to fragmentation (ogi ^ 1). 
However, as shown in the next section, stellar irradiation 
provides a temperature floor that stabilizes self gravity as 
the disk mass decreases. In fact, the self-similar solutions 
above rarely have time to obtain before the disk becomes 
irradiation-supported. 



0.43 Myr *4«:ioA/-Joo(A^d/0.1M* 



(7) 

A reasonable definition for the end of the embedded 
phase and the onset of "isolated" accretion is when the disk 
mass has decreased sufficiently that the accretion timescale 
of the disk tvisc oc M^^ (eq. [7]) exceeds the infall accretion 
timescale from the p rogenitor molecular core tc ~ 10'' years 
(|McKee fc Tanll2003l ). By equating fvisc,gi = tc, we thus find 
that the disk mass at the end of the embedded phase is: 

Md 



0.13Qo*/^yl2(tc/10" yr) 



l/6p-l/6 
100^0,100- 



Robustly then, the disk mass is ~ 10 — 20 per cent of the 
stellar mass at the beginning of the disk's isolated evolution. 

Once the disk no longer accretes appreciable external 
mass and angular momentum, its subsequent evolution oc- 
curs with an approximately constant total angular momen- 
tum Jd oc Md{GMtRdY^'^ , neglecting the effects of disk out- 
flows. As matter accretes, Rd thus inc reases oc MJ^ , i.e. the 
disk viscously spreads to larger radii (|Pringlelll98ll ). 

For times much greater than the initial viscous time 
tvisc,o ~ 10^ yrs, the properties of the disk at Rd asymptote 
to a self-similar evolution which is characterized by f ~ fvisc, 
where tviac is the viscous time at Rd (see 331) • For example, 
from equation we infer that tvisc oc M~[^R^^ oc MJ"*, im- 
plying that Md oc t'^l^"^ ~ As discussed in Appendix 
A, since gravito-turbulent disks obey a non-linear diffusion 
equation, it is not clear that their the viscous evolution is 
indeed controlled by the disk properties at the peak radius 
Rd- Nevertheless, we argue that a self-similar approach may 
still be justified, which results in the following solution (Ap- 
pendix B) 



R-l/6 
*, 100 -"-0,100 



-1/4 



Rd 



2.2Rd,o{t/Uisc,o)^^^ 



T\r^ ~ 130KQo ^''^Atjif Aff''^ 



-4/3 



100-^0,100 



{t/t 



visc.O / 



\rj ~ 230 gem 



H/R\ 



(9) 
(10) 

(11) 
'^/*(12) 

{t/Uis.fi)-"'' (13) 
QmQ-'''n]^,'M:li,lR-\il{t/U,^^,,r^'^ (14) 



-,-4/3 1/6 
CO '^10 



2300Qo"'/'kI^'M^'''^ 



100 

13/6, 



100-^0,100 



■0,100 (^/^visc,( 
5/4 



2.3 Irradiation-Supported Phase 

Irradiation from the star during the embedded stage acts 
to stabi lize proto-stellar disks from gravitational instabil- 
ity (e.g. Matz ner &: Levinll2005l ). As we now discuss, it also 
supplies longer term support, effectively shutting off gravito- 
turbulence soon into the disk's isolated evolution. 

High mass stars radiate at a significant fraction r^cdd of 
the Eddington limit, with a luminosity = r/cddicdd = 
1.3 X 10""'77oddA'/*,ioo ergs s~^0 Stars of mass A/* = 
25(100) AfQ, for instanc e, have 7?cdd = 0.1(0.4) up on entering 
the main sequence fe.g. iMaeder fc Mevnetlll987l ). Assuming 
that the disk is flared, such that the scaleheight H increases 
^ 'with radius faster than oc _R, the stellar flux incident on the 



disk normal at radius R is: 



Lf 
47r7?2' 



(17) 



where / — Tjf{H/R) is the fraction of the total flux inter- 
cepted and absorbed by the disk and rjt ~ few characterizes 
the precise flaring geometry an d number of vertical scale - 
heights to the disk photosphere (jChiang fc Goldreich|[l997l ). 
which depends on details such as grain settling. In our an- 
alytic estimates below we combine r^cdd and rjt into a single 
constant rj = rjcddVi ~ 1- of our results depend sensi- 
tively on Tj. 

Irradiation dominates over viscous dissipation in heat- 
ing the midplane when Fi^ > FvIsc/t, where T^visc is the 
heat flux due, i n this case, to the nominal level of gravito- 
turbulence (e.g. lRafikovll20()9l l. When this condition is sat- 
isfled the resulting midplane temperature Tin is determined 
by equating Fj^r with the vertical radiative flux F ~ SaT'^/S 
(jKratter et al.ll2010h . giving 



686Kn'^'Ml^,l,R-%l 



(18) 



corresponding to a disk scaleheight 



,100 -"-d, 100 



(19) 



"r 

The condition that the disk is irradiation-supported and 
that equations (|18|l and p9p are valid, namely that Firr > 
Fvisc/T, is equivalent to the requirement that Ti^ exceed Tgi 



^ Accretion luminosity is generally negligible in comparison to 
intrinsic stellar luminosity for massive stars. 
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(eq. [3])- This condition is satisfied for disk; masses below 
tfie critical mass 



Ma 



O.I2O0 M277 



1/7, 



-3/7, 
,100-' 



.2/7 
'd.lOO 



(20) 



Note that this mass is comparable to the initial disk mass 
estimated in equation ([8]). This shows that the isolated 
gravito-turbulent stage discussed in ^2.2\ is. at best, short- 
lived because Tgi decreases rapidly (xt~^ during the gravito- 
turbulent phase (eq. [H]) and irradiation provides a stable 
temperature floor. In fact, in some cases Md in equation (|20[) 
exceeds that given in equation in which case the disk 
is irradiation-supported from the very onset of its isolated 
evolution. 

Once the disk midplane becomes supported by exter- 
nal radiation, gravito-turbulent heating is no longer neces- 
sary to maintain Q > Qo ~ 1- Without Gl-induced torques, 
the most viable alternative source of disk viscosity is turbu- 
lence generated by the magneto-rotational instability (MRI; 
iBalbus fc Hawle\Ji 19921 ). If the MRI operates throughout the 
entire disk, the resulting turbulent shear stress can, for our 
purposes, be app roximately described by a Shakura-Sunyae v 
a- viscosity (e.g. iHawlev et"al] Il995l : iFromang et al.ll2007l ). 



with typical values a 



10" 



10 ^ inferred from simu- 



lations and obse rvations (e.g. IChiang fc Murrav-Clavll2007l : 
iKing et al.|[2007l ). The accretion timescale in an irradiation- 
supported a-disk is thus given by 



aCsH 
0.44 Myr(Q/10' 



2/7^5/14 pl3/14 



R 

,iao-"-ii,ioo ' 



(21) 



where we have used H/R from equation (|19|) . Note that 
for values a > 10~^ which are expected for fully-ionized 
disks, the accretion timescale is rather short compared to 
the lifetimes ~ 3 — 10 Myr of very massive stars. 

Proto-stellar disks are, however, dense, cold, and hence 
weakly ionized. When the conductivity is sufficiently low, 
MRI modes are stabilized because the gas collisionally de- 
couples from the magnetic field and ionized species decou- 
ple from the neutral molecular hy drogen that comprises 
the bulk of the disk's mass fe g. iBlaes fc Balbud 1 19941 : 
iReves-Ruiz fc Stepinski 199l|; |jinlll996l ). In particular, the 
surface column through disk near Rd at the beginning of its 
irradiation-supported evolution is: 



1.7 X lO'^gcm" 



-1^1/7,^4/7 



12/7 



!o V ' 100 -^d, 100 



(22) 



where we have used equation (|20|) . This large column shields 
the midplane from ionizing radiation, perhaps creating a 
"dead zone" near the rnidplanc in which little MRI turbu- 
lence is generated ()Gammielll996l : iTurner et al.ll2007l ). 

Determining the presence and extent of the dead zone 
requires estimating the surface column Ea that is sufficiently 
ionized to become MRI active. Cosmic rays are usually as- 
sumed to be the chief ionizing agent at large radii in proto- 
stellar disks becau se they penetrate to a signiii cant depth 
E ~ 10^ g cm"^ jUmebavashi fc Nakanol ligsil ). However, 
the regions around very massive stars may be shielded from 
interstellar cosmic rays by magnetized winds from the star or 
disk. This is particularly true because massive stellar winds 



are more powerful throughout the entire stellar lifetime than 
the outflows from low mass stars during their pre-main se- 
quence phases. 

Assuming that a fraction /b <S 1 of the total luminosity 
of a massive star's wind Lw — (l/2)Mw^'w is carried in the 
form of Poynting luminosity Lp = {B^c/An) x {AtvR'^), then 
the ratio of the Larmor radius of a cosmic ray proton i?L = 
Eci/eB to the radius R is independent of radius: 



^ ^ 10-^ (J^) 
Ed 



-1/2 



( 



GeV 



10-6AfQyr-i 



'1/2 



(l03kms-i) ' ^^^^ 



where Ea is the cosmic ray energy, B is the magnetic 
field, and Uw and iWw ~ 4:Trp„v„R^ are the wind veloc- 
ity and mass-loss rate, respectively. We have scaled Mw ~ 
10"'' - 1O"^M0 yr-i and Vw lO^ '^ - lO^ '^ km s'l to typi- 
cal values f or main- sequence and post-main-sequence O star 
winds (e.g. Ivan Burenlll985l ). 

Equation (|23|l illustrates that cosmic rays near the peak 
of the energy spectrum {Ea ~ 10 MeV - GeV) have Ri^ <C R, 
even for very small values of /b. Interstellar cosmic rays 
should thus act as a fluid and remain exterior to the ter- 
mination shock of the wind with the ISM, which is ^ 10^ 
AU for typical values of the Lw and the ISM pressure. Al- 
though in principle a modest fraction of cosmic rays (espe- 
cially those with higher energies) could diffuse upstream to 
ionize the disk, the fraction that reaches the disk at radii 
R <J 10^ — 10'^ AU is probably small. Cosmic rays could also 
be accelerated locally near the disk, due to e.g. shock accel- 
eration in the stellar wind or disk corona, but the likelihood 
of this possibility and its relative importance are difficult to 
assess. 

X-rays are another source of ionization in pro to-stellar 
disks (|Glassgold et al.lll997l : irgea fc Glassgoidlll999l). star s 
are strong sources of soft X-rays ( Berghoefer et ai1ll996l ). 
with ty pical lumino sities Lx ~ 10'^'^ — 10^^ ergs s~^ (~ 
10~^Z/*; iFlaccomio et al..,2003< ), which are a factor ~ lO'' — 
lO'' larger than the typical X-ray luminosities of T Tauri 
stars. However, the X-ray flux at R ^ 10^ AU is similar to 
that experienced by a T- Ta uri disk at a _R ^ f ew AU . The 
latter case was analyzed by 'igea fc Glassgoldl l|l999l ). who 
found Ea ~ 10 g cm~^, relatively independent of details 
such as the precise X-ray spectrum. Absent cosmic rays, 
we thus take Ea 10 g cm~^ as a fiducial estimate of 
the active layer in massive proto-stellar disks. We acknowl- 
edge, however, that in general the active depth (and even 
its precise definition) depends on a variety of complex and 
interrelated processes, in cluding the precise c onditions for 
the growth of the MRI (|Fleming et all I2OO0I : iPessah et all 
2QQX); the abundance of gas phase m etals; the size, evolu- 
tion, and settling of d ust grains (e.g. iFromang etlo] I2OO2I : 
iBai fc Goodmanll2009l); and th e presence of turbulent mix- 



ing fe.g. lllgner fc Nelson|[200^ . l2008l : [Turner et al.|[2007l ). 

Because the initial column Eq (eq. [22]) greatly exceeds 
Ea ~ 10 g cm~^, proto-stellar disks around massive stars 
likely possess an extensive dead zone at the onset of the 
irradiation-supported phase. Although there have been sev- 
eral studies of [and speculations regarding] the effects of a 
dead zone on disk angular momentum transport, the present 
theoretical picture remains largely incomplete (Sto ne et al.l 
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|2000| ). Below we compare and contrast three qualitatively 
distinct possibilities for the viscous evolution of the disk 
during the irradiation-supported phase. This discussion mo- 
tivates the strength and form of the viscosity that we adopt 
in our disk models in JJS] 



2.3.1 Marginal Gravitational Instability 

One possibility is that angular momentum transport is 
largely ineffective, such that the accretion timescale greatly 
exceeds the stellar lifetime ~ 3 — 10 Myr. In terms of an 
effective a— viscosity (eq. [5T]) this requires a <^ 10~^. Un- 
der this circumstance, one might expect the disk would not 
evolve significantly (absent external influences) after it be- 
comes irradiation-supported. The final mass at stellar core 
collapse would thus simply be given by equation (f20)) . 

This analysis neglects, however, the effects of stellar 
mass loss. Since the Toomre parameter Q is proportional to 
the stellar mass (eq. PP), a decreasing value of M*(f) could 
drive the disk to again become gravitationally unstable at 
later stages of stellar evolution. Because angular momentum 
transport due to gravitational instabilities is rapid (eq. [7]), 
the onset of GI would rapidly deplete the disk's mass until it 
again becomes irradiation-supported. Stellar mass loss thus 
locks a disk that otherwise does not evolve into a state of 
marginal gravitational instability {SQ = 0), such that the 
accretion rate is tied to the stellar mass loss rate by the 
relationship 

MdlsQ^o = (2Afd/llM,)M*, (24) 

where we have used equation H20p and assumed that the 
disk spreads with constant total angular momentum Jd oc 
MdiGM^^Rdf^. Thus, in this case the final mass at stellar 
core collapse is smaller than its initial value at the onset 
of irradiation-support by a factor (M^j/M-kfi)'^^^^, where 
M*j = Mt,{t = tiifc) is the final stellar mass. The final ra- 
dius is Rdj/Rdfi = {M^,J /Mt,fi)~^^^^^ . Unsurprisingly, neg- 
ligible non-gravitational viscosity results in the most mas- 
sive, compact disks which can survive until stellar death. 



2.3.2 Active Zone Draining 

Another possibility is that accretion does occur in 
irradiation-supported disks, but it is restricted to regions 
of th e disk which are sufficiently ionized to be MR I active 
(e.g. iGammid [l99^ : IChiang fc Murrav-Clavl [20071 ). In ad- 
dition to the externally-ionized surface layer Sa previously 
discussed, the entire disk is MRI active due to the coUisional 
ionization of alkali m etals when the midplane temperatur e 
exceeds Ta « 10^ K l|Gammiel [l99^ : iFromang et al.ll2002h . 
Using equation (|18p this condition is satisfied interior to the 
radius 

Ra = R{T = To) ~ 40 AU (r„/10^ K)-'"\'^'''m]%^, 

(25) 

From equation (|21|) . the accretion timescale at Ra is thus 
fviBc.irrlfl^ ~ 0.02 Myr {a/Q.l)S'''Hll%o{Ta/lG' K)-"''^ 



where we now scale the viscosity to a larger value a --^ 0.1 
which is appropriate for the MRI-active inner rim of the 
disk. 

Although tvisc,irr|js„ is vcry short, it is not the timescale 
for the entire disk to accrete. At radii R > Ra only the ac- 
tive surface column Ea feeds the inner disk and, as a re- 
sult, the a ccretion rate reaches a minimum just outside Ra 
l|Gammie|[l996 l. The inner edge of the dead zone thus acts 
like a spigot by controlling the net inflow rate for the entire 
disk. This results in an effective accretion timescale for the 
disk to "drain" through its active layers tdrain which is larger 
than equation H26p by the ratio of the total disk mass to the 
"active" mass at i? « Ra, Ma ~ 7ri?^Ea: 

tdrain ~ {Md / Ma)t 

visc,irr | Ra 

^ 32Myr(a/0.1)~\~^ x 
M.,ioo(Md/0.1M*)(T,/10=* K)'/'(Sa/10gcm-")-\ 

(27) 

For typical values of the disk's mass at the onset of the ir- 
radiated phase (eq. [20]), tdrain generally exceeds the stellar 
lifetime tufc ~ 3 - 10 Myr for M*,o > 25Mq. This suggests 
that the disk is unlikely to evolve significantly prior to stellar 
death if it accretes solely through its active regions; its evo- 
lution would thus revert to the case of negligible irradiation- 
supported viscosity described in i)2.3.1l In the next section 
we discuss the perhaps more plausible situation that limited 
turbulent angular momentum transport (perhaps seeded in 
the active layers) acts throughout the entire disk. 

2.3.3 (Reduced) Global Alpha Viscosity 

A final possibility is that proto-stellar disks undergo a vis- 
cous evolution that is qualitatively similar to that of fully- 
ionized disks, but with a lower "effective" value of the a pa- 
rameter (yet larger than the case of negligible viscous evolu- 
tion described in i]2.3.ip . This could occur if MRI turbulence 
is limited to the active s urface layer Eg, but thi s still gener- 
ates a Umited R eynolds ( iFleming fc Ston ^120031 ) or Maxwell 
st ress ([Turner et al. 2003) midplane. The simulations 

of [Fleming fcStond (|2003[ ). for instance, found that the ef- 
fective angular momentum transport was reduced by a fac- 
tor of ~ 10 in the case of an active layer column which was 
~ 20% of the total column. In this case we would predict 
"effective" viscosities a ~ 10~^ — lO"'^ assuming "standard" 
(fully-ionized) values a ~ lO"'^ — 10~^, at least when the 
disk surface density obeys E|i{^ 2> Ea . In general we might 
expect a to be a decreasing function of E/Ea (which asymp- 
totes to the fully ionized value of a as E — >■ 2Ea), but since 
this dependence has yet to be determined we hereafter as- 
sume a — constant. 

Assuming a constant a disk model with an initial mass 
as given in equation (|20p . the late-time similarity solutions 
(analogous to eqs. [9] — [15]; see Appendix B) are given by 

Md ~ 14M0Qo S'^'<''ioo^o:?oo(iAvisc,o)-'/" (28) 

Rd ~ 1.7i?d,o(t/tviBc,o)"''" (29) 

T\r, ~ 550Kr;^/^A4%i?-?/o'(t/t.iec,o)-'/" (30) 

E|b, ~ 200gcm-"Oo"''^ioA<''ioo^o"wo'(iAvisc,o)"''''" (31) 

t\r^ ~ 2OOOQ0 'KioM,*''/ooi?-i^/^(t/tvi.c,o)"''''" (32) 



(26 
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H/R\a, ^ 0.0777^/^Af-i^/oi?o;ioo(iAvi=c,o)*/" (33) 

OIb, =i 0.5Qo(tAvisc,o)"''", (34) 

where fvisc.o here refers to equation (|21|l evaluated at 
the beginning of the irradiation-supported phase, and we 
have taken A — Ai^ ~ 6.8, as appropriate for an 
irradiation-supported a— disk in steady-state (see Appendix 
A). These solut ions are analogous to those obtained by 
iHartmann et al. I 11998) as applied to T Tauri disks. 

From the above note the following: (1) the relatively 
slow evolution of the disk oc t~''^^^ implies that the disk 
can have substantial mass on timescales S> tviac.o; (2) our 
assumptions that the midplane is optically thick and that 
the temperature remains sufficiently high for our assumed 
opacity (^ 100 K) remains valid for several viscous times; 
(3) Q increases with time, implying that the disk evolves to 
an increasingly gravitationally-stable state (i.e. away from 
the marginally unstable state described in i|2.3.H) : (4) the 
disk column E remains > Ea ~ 10 g cm^^ for at least a 
few viscous times, thereby justifying our use of a "reduced" 
value for the viscosity a on timescales J> tvisc.o. 



3 DISK EVOLUTION MODELS 

We now combine the results and intuition gained from the 
previous section to construct a simple model for the iso- 
lated evolution of massive proto-stellar disks following the 
embedded stage. Although several detailed time-dependent 
models of proto- stellar disks have been constructed by previ 



models or proto- stellar disks nave been constructed by previ- 
ous authors (e.g. Clarke et al.ll200ll:lMa,tsuvama et a,l.ll2003l : 



lAlexander et all bood : IZhu et all l2009l : ICorti et all l2009f) . 

our calculations are unique in focusing on the evolution of 
disks around very massive stars for the entirety of the stellar 
lifetime. The time-dependent disk properties we obtain are 
used to assess the effects of external, destructive influences 
in El 

3.1 Self-Similar Ring Model 

Our d isk model closely follows that used bv lMetzger et all 
(|2008l ') t o study accretion follo win g compact object merg ers 
(see e.g. ICannizzo et al.|[l990l and IHartmann et al.lll99a for 

related self-similar models). 

At a given time t, the disk can be divided into three 
regions, depending on the local viscous time tvisc, which 
generally increases with radius for a fixed total angular mo- 
mentum (e.g. eq. [5T]). At small radii, tvisc < t, and the disk 
enters a steady state with M oc i/T, constant, where i/ is the 
kinematic viscosity. Larger radii where tvisc ~ t contain the 
majority of the disk's mass and angular momentum. This 
region determines the viscous evolution of the rest of the 
disk, including the mass accretion rate that is fed into the 
interior, steady-state region. Exterior to this point is a re- 
gion with tvisc > t, but this contains only a small fraction of 
the mass and does not significantly affect the viscous evolu- 
tion (see, however, the discussion of gravito-turbulent disks 
in Appendix A). 

Our model treats the disk as a single annulus (the 
"ring") that is evolved forward in time. The properties of 
the ring, such as the surface density E and temperature T, 
are representative of the location R ~ Rd where the local 



mass Er peaks. The time evolution of the disk is deter- 
mined by the conservation equations for mass and angular 
momentum: 



d 



- (AttEJ?!) = -Md, 



- [B{GM.Rdy^\T.Rj] = - j = 0, 



(35) 
(36) 



where Md is the mass accretion rate, J is the angular mo- 
mentum loss rate, which we take to be zero (i.e. we neglect 
disk outflows), and A and B are constants which relate the 
total disk mass and angular momentum to their local values 
near Rd- At early times the constants A and B depend on 
how matter is initially spatially distributed. At times much 
greater than the initial viscous time, material spread in a 
manner determined by the viscosity. As described in Ap- 
pendix A, A and B are determined by setting the solution of 
our simplified ring model at late times equal to the solution 
for a spreading ring. In the irradiation-supported accretion 
phase f ^2.3[) the viscosity obeys u oc r^^'^^* and this analy- 
sis results in exact values (Aim Birr) = (6.80,5.85). In the 
gravito-turbulence phase ( t|2.2|l the situation is more compli- 
cated; since the viscosity obeys v oc E^r^^, the surface den- 
sity evolves via a non-linear diffusion equation which con- 
centrates the mass approximately equally per unit decade in 
radius (see Fig. lAip . As described in Appendix A, we take 
Agi ~ 2 in this case. 

The accretion rate depends on the characteristic mass 
and viscous timescale of the ring as 



Md = fMd/Uu 



(37) 



where the factor / is set, like A and B, to match the solution 
of a spreading ring (Appendix A) and has the values /gi ~ 3 
and fin ~ 1.4 for the gravito-turbulent and irradiation- 
supported phases, respectively. The viscous time tvisc is 
determined using equation ([7} when the disk is gravito- 
turbulent (i.e. when Tgi > Tirr; eq. [21]). When the disk is 
irradiation-supported (i.e. Tirr > Tgi) we use equation (|21|l 
assuming a fixed value a ~ 10^^ — lO"'^, as motivated by 
the discussions in i]2.3.3l and ^2.3.21 

Equations (|35p and (|36p provide two coupled equations 
to be solved for the dependent variables Rd and Add ~ 
AtvRIT, given i?d(0), Md{0) and M^,{t). The initial condi- 
tions reflect those at the end of the embedded phase. In 
particular, we take the initial disk mass as the minimum 
of equations (|8]) and (I20|l . assuming a value for the critical 
Toomre stability parameter Qo = 1. In all of our models 
we take an initial radius Rd,o = 100 AU, as motivated by 
the work of KM06 and the discussion in f\2.1\ We smoothly 
interpolate the parameters y — A, B, and / between the 
gravito-turbulent (gi) and irradiation-supported (irr) accre- 
tion regimes by employing the simple prescription 



y = j/irr exp[-rgi/Tirr] + 2/gi(l " exp[-rgi/Tir 



(38) 



For all stellar masses we take tc = 10^ years for the duration 
of the embedded phase, although our results are insensitive 
to this choice. We run the calculation until the end of the 
stellar lifetime {t — tiifc). 

Stellar evolution models show that wind mass- 
loss is important throughout the main sequence and 
post main sequence evolution of very massive stars 
(e.g. iMaeder fc Mevnetu 1987 1. Our models take into account 
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the efTects of a decreasing stellar mass Mi,(t), both through 
its affect on disk properties (e.g. fvisc) and in the angular mo- 
mentum equation fe g. 1361). In most of our m odels we use 
M*(t) and from lMaeder fc Mevnetl (| 19871 ). but we also 
perform comparision calculations in which we neglect mass 
loss entirely, as may be justified for stars with low metallicity 
or if the mass loss rates employed in standard stellar evolu- 
tionary models are overestimated due to e.g. wind clumping. 
However, for simplicity we neglect the resulting shortening 
of the stellar lifetime. 



3.2 Results 

We calculate evolutionary disk models for stars with ini- 
tial (post-embedded) masses A/* = 25 — 12OM0 and for 
different values of the viscosity a during the irradiation- 
supported phase. Table [1] summarizes the properties of the 
disk at stellar core collapse [t — tufc). For the case of very 
low a 10~^), the disk evolves in a state of marginal grav- 
itational instability, as described in i]2.3.1l As our results in 
Table [1] confirm, this results in the most massive, compact 
disks at core collapse. 

Sample solutions for the case of finite values of a are 
shown in Figures [T] and [2] As we now describe, these cases 
highlight important differences between the evolution of 
disks around less massive and more massive stars. 



3.2.1 M*,o = 25Mq 

Figure [T] shows the disk properties at the peak radius Rd as 
a function of time for a star with initial mass i\f*,o = 25Mq 
for an assumed viscosity a — 10~^ during the irradiation- 
supported phase. As expected, the disk experiences a very 
brief gravitationally unstable phase of rapid evolution ( i)2.2p . 
after which it becomes irradiation-supported and its viscous 
evolution slows considerably ( H2.3p . The disk mass at the 
beginning of the isolated phase is ~ 5Mq and this decreases 
to ~ 2Mq by core collapse at t = tufc = 8.3 Myr. As a 
result of angular momentum conservation, the disk spreads 
to 1200 AU hy t = tiifc, with a significant fraction of 
this expansion occurring in just the last ~ 10^ years, when 
the stellar mass decreases to ~ UMq due to stellar winds. 
This late spreading occurs even absent significant accretion 
because conservation of total angular momentum requires 
that Rd increase oc M^^ as M* decreases, even at fixed disk 
mass. Indeed, in the otherwise identical model neglecting 
stellar mass loss, the final radius is only ~ 700 AU (see 
Table [g. 

We also note that the surface density E decreases below 
10 g cm~'^ for t >4 Myr. Since this is similar to the surface 
layer Ea which is MRI active due to X-ray ionization ( i]2.3|) . 
the dead zone in the disk midplane may be eliminated. A 
larger value of a than the fiducial value 10~^ may thus be 
more physical at these late times. A model calculated with 
a = 10~^, however, results in a much larger final radius ~ 
10* AU (Tabled . As we discuss in 331 such an extended disk 
is particularly susceptible to photo-ionization mass-loss and 
stellar collisions. This illustrates that massive disks may be 
unlikely to survive ar ound 25M(^ stars, a result in agreement 
with other work (e.g. lGorti et al.ll2009l ). 
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Figure 1. Time evolution of the disk properties for a star with 
initial mass M*^o = 25Mq, calculated assuming a value a = 10~^ 
for the viscosity during the irradiation-supported phase. The disk 
properties shown include the total mass and the peak radius 
Rd; and the surface density E, midplane temperature T, Toomrc 
parameter Q, and the scaleheight H/R, all evaluated at the radius 
Rd ■ Also show n is the stellar mass Mj, (t) from the stellar evolution 
calculations of lMaeder fc Mevne"3 |l983). 
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Figure 2. Same as Figure[T] except calculated for an initial stellar 
mass M*,o = 85Mq. 



3.2.2 M*,o = 85M0 

Figure [2] shows the disk properties for an initial stellar mass 
M^, — 85Mq and a — 10""^. In this case the disk mass at the 
beginning of the isolated phase is ~ 17AfQ, which decreases 
to ~ 8Mq by core collapse. Importantly, the surface density 
remains > 10 g cm~^ ^ Ea throughout almost the entire 
stellar lifetime (as opposed to the Af* = 25AfQ case), despite 
significant spreading due to stellar mass loss. A dead zone 
may thus be present until stellar core collapse, implying that 
the "reduced" value of the viscosity that we have adopted 
is self-consistent. 

Given that more massive stars have shorter lifetimes 
and longer-lived dead zones, we conclude that they likely 
possess more massive, compact disks at core collapse when 
considering isolated accretion alone. 
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Table 1. Properties of Relic Proto-Stellar Disks at Core Collapse 





life 








Rd 






(*f0) 


(Myr) 








(AU) 




(g cm-2) 


120' 


3.5 


~ 


64.0 


13.3 


240 


0.09 


330 


120 


- 


10-=* 


- 


10.4 


670 


0.14 


12 


120* 


- 


10"' 


- 


10.2 


380 


0.09 


35 


120 


- 


10-2 


- 


4.4 


3800 


0.22 


0.13 


120* 


- 


10-2 


- 


4.2 


2300 


0.15 


0.34 


85t 


3.8 


~ 


38.4 


9.4 


300 


0.12 


150 


86 


- 


10-' 


- 


7.7 


830 


0.18 


5.8 


86* 


- 


10-" 


- 


7.4 


400 


0.10 


22 


86 


- 


10-2 


- 


3.1 


5200 


0.30 


0.05 


86* 




10-2 




2.9 


2700 


0.17 


0.17 


60+ 


5.4 





21.3 


7.4 


410 


0.18 


62 


60 




10-3 




5.4 


1140 


0.24 


2.2 


60* 


5.4 


10-3 




5.3 


440 


0.12 


13 


60 




10-2 




2.1 


8000 


0.42 


0.014 


40t 


5.4 


~ 


9.8 


5.5 


680 


0.29 


17 


40 




10-3 




3.6 


1900 


0.37 


0.5 


40' 


5.4 


10-3 




3.5 


500 


0.14 


6.0 


40 




10-2 




1.3 


13000 


0.67 


0.002 


25t 


8.3 


~ 


14.0 


4.9 


220 


0.17 


140 


26 




10-3 




2.1 


1200 


0.27 


0.64 


26* 




10-3 




2.0 


700 


0.18 


1.8 


26 




10-2 




0.70 


11000 


0.51 


0.002 



Initial (post-embedded) stellar mass; Stellar lifetime; 

Viscosity during the irradiation-supported phase; Stel- 
lar mass at core colla pse (from the evolutionary models of 
iMaeder fc Mevnet] [ToSTh : + Calculated assuming negligible angu- 
lar momentum transport during the irradiation-supported phase 
(see ^2.3.11 1 'Calculated assuming negligible stellar mass loss. 



4 DISK DISPERSAL PROCESSES 

Our results in ij3]show that isolated evolution generally re- 
sults in a relatively massive proto-stellar disk at core col- 
lapse. However, it is well known that accretion alone can- 
not explain the observed lifetimes of disks around lower 
mass proto-stars. In this section we thus address external 
processes that may act to disperse the disk. These include 
photo-evaporation ( i]4.ip . stellar collisions ( i|4.2[) . stripping 
by continuous stellar winds (i )4.3p . and explosive (LBV- like) 
stellar eruptions ( i|4.4[) . In i|4.5l we combine our conclusions 
to assess whether and under what conditions a proto-stellar 
disk is most likely to survive until core collapse. 



4.1 Photo-evaporation 

Perhaps the most important means of disk dis- 

persal is photo -e vapor at ion due to UV irradiation 

llHoUenbach et al.l Il994 ': 'Shu et al.1 1 19931 : I Clarke et all 



I2OOII : iGorti fc HoUenba ch 2009). Very massive stars have 
large Lyman continuum photon luminosities of ~ ]^q48-49 
s~^, which ionize and heat the upper surface layers of the 
disk to a temperature ^ 10* K, somewhat analogous to 
an H II region. This heating drives a Parker-like thermal 
wind from the disk's surface. A critical radius for photo- 
evaporation Rg occurs where the escape speed of the disk 
equals the sound speed of 
= Cs E ~ 10 km s~^: 



10'' K gas, Cs(r = lO^'K) 



i?g = {2GM^/clg) ~ 1800AUAf*,ioo (39) 

iHollenbach et all l| 19941 ) show that the total photo- 
evaporation rate is dominated by outflows from radii near 
J?g. For disks that extend to R > Rg they find that the total 



mass loss rate i^ 

Mph ~ 4 X lO-^M0yr-i(<E>i/lO*''s-i)^/2A4%, 
resulting in an evaporation timescale tph given by 
fph = Md/Mph ~ 

3 X 10'^yr(Md/0.1M,)Af'/oo($i/10''''s-^)" 



1/2 



(40) 



(41) 



Since is generally <^ the stellar lifetime tufc ~ 3 — 10 
Myr, this implies that if a disk spreads to radii ^ -Rg, it will 
photo-evaporate prior to core collapse. 



highly suppressed (Afph oc exp[—R/2Rg], approximately) 
because the sonic radius ~ Rg occurs far out of the mid- 
plane, many sc ale-heights down t he disk's exponential at- 
mosphere (e.g. I Adams et al.|[20o3 ). Disks with outer radii 
Rd less than a few times Rg thus have tph > info and should 
not photo-evaporate. 

For Af* = 25Mq the photo-evaporation radius Rg is 
only ~ 500 AU. This is less than or comparable to the size of 
the disk at core collapse that we found for Af* — 25Mq in 33] 
(Fig. HI Table[T]), even in the models that neglect stellar mass 
loss. Photo-evaporation may thus be effective at dispersing 
the proto-stellar disks of 25Af0 stars. 

By contrast, for Af* = 85AfQ the critical radius Rg « 
1500 AU is significantly larger than the disk radius at core 
collapse for a < 10"^ (Fig.0. While it is true that Rg oc Af* 
may decrease significantly due to stellar mass loss, most of 
this occurs in just the last ~ 10"" years before stellar death. 
Although the photo-evaporative mass-loss rate from the disk 
may thus become large (e.g. ^ IO'^A^q yr"'^) prior to core 
collapse, there is insufficient time to photo-evaporate the 
disk. 

To summarize, although photo-evaporation may dis- 
perse protostellar disks around less massive stars, higher 
mass stars become increasingly difficult to photo-evaporate 
because their deep gravitational potentials imply large val- 
ues of Rg relative to the final disk radius Rd- This statement 
holds true best for disks with small initial radii (as form from 
low angular momentum cores) and for lower stellar mass-loss 
rates, as would preferentially occur for low metallicity. Nev- 
ertheless, we emphasize that significant photo-evaporative 
mass-loss from the disk likely does occur and this may have 
observable consequences for the supernova (see ^6.2^ and the 
appearance of its progenitor ( ij6.ip . 



4.2 Stellar Encounters 

Massive stars generally form in dense stellar clusters (e.g. 
Ide Wit et al]|2005l : iTanI 120071 '). Close stellar encounters are 
thus common and may periodica lly strip off the outer 
edge of the proto-stellar disk (e.g. IClarke fc Pringlel Il993l : 
iHellej [l995l : iHaU et al.l 119961 ). Although the most massive 
clusters in the Milky Way have central stellar densities 



Since the outflows from ma s sive s tars obey the "strong wind" 
condition of iHoUenbach et al. I 1I1994I) . the mass loss rate in equa- 
tion 1 140 I I formally applies only to disks with outer radii Rd ~ Rg ■ 
In the present situation this distinction is, however, not essen- 
tial; because the evaporation timescale tpji (eq. 1411 ) is generally 
shorter than viscous spreading time, photo-evaporation will trun- 
cate the outer edge of any disk that expands to radii > Rg. 
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rii, ~ lO"* — 10^ pc ^ (e.g. iMcCaughrean fc Stauffeij[l993 : 
iGarmire et al.|[2000l ). the density around a typical O star 
througho ut its lifetime d epends on uncertain detail s such as 
the IMF jAdams fc My ers 2001; Massi et aLll2006l) and the 
degr ee of primordial (e.g. .Huff fc Stahlerii2006l ) and dynam- 
ical (jPortegies Zwart et al.ll2004l ) cluster mass segregation. 

In the limit of strong gravitational focusing, the colli- 
sion timescale for stars with density n* and one-dimensional 
vel ocity dispersion a with a di sk of cross section tvR'^ is given 
by (|Binnev fc Tremaine|[l987l ): 



4.3 Stellar Wind Stripping 

Mass loss due to interaction with the stellar wind is an- 
other wa y to dispe rse the disk telmeg reenlflOTSl : lYorkell2004l : 
iMatsuvama et al.l [2009i ) . I n evaluating this possibi lity we 
closely follow the results of IMatsuvama et al.l (|2009l ). here- 
after M09, who consider a model in which the central wind 
obliquely strikes the flared disk surface, driving material out 
of the system by entraining it in a mixing layer along the 
disk surface. M09 estimate that the timescale for significant 
mass-l oss from the disk due to stellar wind stripping is given 
= 4xl0^yrK/10*pc-^)-'7?d,looW3kms"')"'A/,;i'oo-(42) by (cL iHollenbach et al.|[2000h 



Although this timescale is very short, most collisions occur 
with stars near the peak of the IMF, which have a typical 
mass M*,o ~ 0.5Mq. Yet, because incoming stars are typ- 
ically on near-paralDolic [zero-energy] orbits (e.g. ,Ostrikeii 



ically 



by energy conservation alone only stars with masses 
^ Md are capable of unbinding a significant fraction of 
the diskQ For most impact trajectories the disk mass re- 
moved w ill be s i gnifi cantly less than that of the colliding 
star (e.g. iHellei] 09951 '). 

We may crudely estimate the timescale tsc for the disk 
to be depleted by stellar collisions if we assume that colli- 
sions with stars of mass M*,c dominate the mass loss and 
that the average mass ejected in a collision is a fraction rjsc 
of Mt_c. We accomplish this by replacing the stellar density 
n* in equation (|42|l with just the fraction (M*,min/M*_c)^ '^ 
of the total density n*,tot in stars with mass ~ Af*,c, 
where Mi-^min ~ 0.5Mq and we assume a Salpeter IMF 
dN^/dM^ (X M-^ '-^ above M*,min. This gives a depletion 
timescale 



tsc — tc 



i 104 pe- 



rl*. tot 



O.IM* 

Rd 



Md 

0.3 



We have normalized rjse. to a value 0.1 which has been conser- 
vativ ely extrapolated from numerical simulations (e.g. lHellerl 
[l99i) to the case in which the disk mass Md and colliding 
stellar mass M*.c are comparable. Simulations suggest that 
?7 decreases moving to larger ratios M£j/Af*,c, but to our 
knowledge no calculations of the disk mass loss have yet 
been performed in the regime Af*.c > Md- 

Equation H43p shows that tse may be comparable to the 
stellar lifetime for disks with masses Md ~ 0.1 Af* and radii 
Rd ~ 100 AU in clusters with < 10* pc"^. Though this 
estimate is extremely crude and the importance of collisions 
depends sensitively on the conditions specific to an indi- 
vidual massive star, it appears plausible that massive disks 
could survive until core collapse in "typical" cluster environ- 
ments. We note, however, the obvious fact that than since 
tae oc J?^^ more extended disks become increasingly suscep- 
tible to dispersal. 



4 For example, in the extreme case that a star of mass Mi, is 
captured into a circular orbit at radius Rd, only a disk mass 
Md = Mi, could become unbound due to the resulting energy 
release. 



Afdcos /3cs(l-h?^ 



AfwWw(^ sin^ 7 

--(^)"(ti^)(o^I:)( 



Mi, 



\ WOMq 



UO^kms 



Af« 



IO-sA/q yr-i 



sm 7 
10-2 



(44) 



where f < 0.1 is the Mach number of the disk material that 

^ i 

mixes into the shear layer ( [Canto fc Raga 1991), Cg is the 

sound speed in the disk below the mixing layer, and P/y are 

the angles that the mixing layer makes with the midplane 

and the wind, respectively (see Fig. 1 of M09 for the relevant 

geometry). In deriving Equation (|44|) . the disk density p used 

in determining the mass-flux into the mixing layer S oc p^Cs 

is obtained by equating the incident wind ram pressure oc 

MwVw sin^ 7-R^^ with the disk thermal pressure oc pc^. 

In the second and third line in equation (|44l) we have 
assumed that ^ <C 1 and cos/3 ~ 1, and we scale sin^ 7 to 
the value ~ lO"'^, motivated by the modest wind- interface 
incidence angles 7 ~ 1 — 5° found by M09 as a part of self- 
consistent numerical calculations (see their Fig. 5). Although 
M09 consider lower mass loss rates and velocities than those 
appropriate for massive stellar winds, they find that 7 is 
relatively insensitive to the wind properties (if anything, 7 
appears to decrease with the wind power; see their Fig. 6). 

The midplane temperature in the disk (T > 100 K) sets 
a lower limit to the sound speed Cs ~ 1 km s~^. Thus, for 
properties typical of O star winds, Af„ ~ 10^^ — W^Mq 
yr"-^ and ~ lO^-^ - 10^'^ km s~^, equation shows 
that Tws generally exceeds the stellar lifetime fufe ~ 3 — 10 
Myr. The M09 model does not consider the possible entrain- 
ment of matter from the inner rim of the disk at 7? ~ 7?a, 
where the wind intercepts the disk head-on (7 ~ 90°) and 
the mass entrainment rate E oc sin^ 7 could in principle be 
much higher. However, the total mass entrainment rate in 
this case A^ws oc EA oc Acs must also be reduced due to the 
smaller surface area A of the rim relative to the disk annulus 
(by a factor ~ H/R\r^ ~ 0.1) and the higher sound speed of 
the inner rim. All things considered, the qualitative conclu- 
sion that Tws ^ Tiifc for massive stars remains intact. More 
detailed calculations, especially a more precise determina- 
tion of the "mixing efficiency" ^, will however be required 
to confirm this conclusion. 
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4.4 Giant LBV Eruptions 

A final means to disperse the disk is via singular, ener- 
getic eruptions from the star, such as those which may oc- 
cur d uring the luminous blue variable (LBV) phase (e.g. 
IContilll9 84: Boha nnanlll997h of stars with masses > 40Mq 
( Langer et al.,, 19941 '). Assuming that the mass of the disk is 
concentrated near its outer edge, the gravitational binding 
energy of the disk is approximately 



2R 



10^^ ergsMl^oo{Md/0.lM^)R, 



'd,100- 

(45) 



Because a fraction ^ H/R\r^ ~ 0.1 of a spherically- 
symmetric outflow intercepts the disk (see Figs. |3]and[Sl), 
an eruption with energy > 10** ergs would in principle be 
sufficient to unbind a disk with mass Aid ~ IOMq from 
a 100 Mq st ar. In particular, 77— Carinae's giant outburst 
in the 1840s (j Smith et al.!!2003). with an estimated energy 
~ 10^" ergs (jSmithl I2008ai ). could destroy even a massive 
proto-stellar disk. 

Whether the disk in fact becomes unbound depends, 
however, on how efficiently the LBV outflow shares its ki- 
netic energy with the disk. To see why the disk-outflow cou- 
pling may in general be low, it is useful to draw an analogy 
with a similar physical system: a supernova s hock impact- 
ing a molecular cloud (iMcKee fc Cowid 1 19751 : iMcKee et al.l 
ll978l : lKlein et al.lll994l ). When the LBV ejecta (with mass 
Mej and expansion speed Vej) reaches the disk, a bow shock 
will form around the disk's inner edge. The resulting pres- 
sure will drive a strong shock through the disk midplane at a 
given by the shock jump conditions 
, where 



1/2 



velocity Vnh ^ Ve\X " 
iMcKee fc Cowielll975l ) 



Pd ^ ( H/R\nA ~' [mA 



(46) 



is the ratio of the midplane density of the disk pd ~ 
Md/2-KRdH\Rd to the density of the LBV ejecta pcj ~ 
M,j/(4^i?i/3). 

If the disk is more dense than the ejecta (x > 1), 
the time that the shock requires to cross the disk radi- 
ally tcr ~ Rd/Vah (the "cr ushiug time," using the termi- 
nology of iKlein et all 1 19941 ) is a factor ^ y^^"^ longer than 
the expansion time the outflow takes to flow around the 
disk toxp ~ Rd/vcj- Thus, when x S> 1 the ram pressure 
of the outflow is not sustained sufficiently long to acceler- 
ate the entire disk to the post-shock speed v^i^, much less to 
the eje cta speed (which re quires several additional crushing 
times; iMcKee et al.|[l987l ). As a result, only a small por- 
tion of the outflow's kinetic energy is imparted to the disk. 
A sufficiently dense disk should thus remain bound to the 
star, even if its binding energy is much less than the kinetic 
energy of the colliding ejecta. 

During its giant eruption, r;— Car ejected a mass Mcj ~ 
IOMq in less than a decade. In this case x ^ t for Md <, 
Mq, implying that e ven a moderately massive disk could 
become unbound f see iDwark ada s fc Balick 1998 ). We note, 
however, that 77— Car may represent an extreme case and its 
eruptions may not b e representative of those from LBVs as 
a whole l|Vinkll2009l ). Observations of the nebulae around 
LBVs typically indicate a wide distrib ution in the ejecta 
mass and energy (Moj ~ 0.1 — IOMq; e.g. lHutsemekers et al.l 




000 



Rd (AU) 



Figure 3. Time evolution of proto-stellar disks in the space of 
disk mass Md (left axis) and radius Rd (bottom axis) for initial 



stellar masses in the range 



25 



12OM0, calculated as- 



suming that the disk resides in a state of marginal gravitational 
stability at Qlfi^ = Qo = 1 ( ^2.3.11 1 throughout the entire stellar 
lifetime, such that the disk spreads solely as the result of stellar 
mass loss. Triangles show the initial location of the disk follow- 
ing the onset of its irradiation-supported phase, assuming that 
the initial radius of the disk is ~ 100 AU independent of stellar 
mass (KM06; see ^2.1|l . Diamonds show the final disk position at 
core collapse. The point at which the radius of the disk equals 
the gravitational radius Rg (eq. [39]) is shown with a dashed line. 
Radii of constant surface density E = Md/A^^TvR^ are shown with 
dotted Unes for S = 10 and S = 100 g cm~^. 



ll994l : ISmith fc Owockilliooj ) . but in most cases it is unclear 
whether th e ejecta primarily results from one or several giant 
eruptions (ISmith fc Owockil 20061) or in a more continuous 
wind (e.g. ICarcia-Segura et all 19961) . As discussed in t)4.3l 
the latter case is unlikely to disperse the disk. Even if LBV 
outflows are dominated by explosive events, we suspect that 
disks with masses > Mq may be sufficiently dense (x 2> 1) 
to remain bound, although hydrodynamical simulations will 
ultimately be required to precisely delineate how large x 
must be for disk survival. 

Incidentally, we note that were an LBV eruption 
to impact a relic proto-stellar disk, the resulting event 
(though not a supernova) would still be relatively luminous 
and could be contributing to the observed population of 
lower-luminosit y transients known as "supernova impostors" 
(|van Dvkll2005l) . 



4.5 Summary: Conditions for Disk Survival 

Taken together, our results suggest that massive disks may 
preferentially survive until core collapse around the most 
massive stars. One reason is that the most threatening dis- 
persal processes, photo-evaporation ( H4.ip and stellar colli- 
sions (i ]4.2|l . become increasingly effective for more radially- 
extended disks. As shown in ^ the disks around more mas- 
sive stars tend to be more massive and compact at core col- 
lapse due to their higher initial masses; shorter stellar life- 
times; and less efficient angular momentum transport dur- 
ing the irradiation-supported phase, which results both from 
their higher surface densities and shielding from cosmic ray 
ionization by powerful stellar winds. Everything else being 
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equal, the more massive disks around massive stars are also 
less susceptible to destruction via stellar encounters and 
stellar wind stripping (simply because more mass requires 
longer to remove) and to giant LBV eruptions due to the in- 
efficient coupling between a dense disk and lower-mass LBV 
outflows. Important caveats to these conclusions include the 
large uncertainties in how the average cluster environment 
and the strength of LBV eruptions depend on stellar mass. 
It is nevertheless plausible that the average disk mass at core 
collapse is a steeply increasing function of stellar mass. This 
is consistent with the observation that long-li ved massive 
disks are rare around B a nd late-type O stars (jNatta et al.l 
I2QOOI : I Fuente et al.l [200^ . but suggests that disk survival 
should become more likely moving to early-type O stars. 

We emphasize that we do not propose that disk sur- 
vival fraction approaches unity moving to the highest mass 
stars; rather, we are pointing to a trend within the confines 
of our assumptions. There are still many ways that even 
disks around very massive stars may be destroyed, one of 
the most severe being the effect that a massive binary com- 
panion would have on the disk, a possibility which we have 
neglected from the onset, but which could nonetheless be a 
crucial factor in the majority of systems. Indeed, if massive 
disks survived around all of even the highest mass stars (e.g. 
Mi,fi > 80Mq), the resulting number of VLSNe would be 
too high. 

We encapsulate some of our conclusions in Figure [3l 
which shows the path taken by proto-stellar disks from the 
beginning of their isolated evolution until core collapse in 
the space of disk radius and disk mass for an initial ra- 
dius i?d,o ~ 100 AU. For concreteness we assume that an- 
gular momentum transport is negligible when the disk is 
irradiation-supported, such that the disks remains in a per- 
petual state of marginal gravitational instability ( H2.3.ip in 
which accretion and disk spreading occurs solely as the re- 
sult of stellar mass loss. For comparison we show the location 
where the disk radius equals the gravitational radius Rg for 
photo-evaporation (eq. [39]) and contours of constant sur- 
face density S = 10, 100 g cm~'^. Figure [3] illustrates that 
disks around the highest massive stars spend the majority 
of their lives at radii <^ Rg and with high surface densities 
(S> Ea ~ 10 — 100 g cm~^), i.e. those conditions which are 
most favorable for survival. 



Pre-Supemova Disk and Winds Xvind (disk)~ 1 0- 1 00 km/s 

* * ,, 4 



5 CIRCUMSTELLAR INTERACTION IN 
CORE COLLAPSE SUPERNOVAE 

Our results from Sections[3]and[4]suggest that when a proto- 
stellar disk survives until core collapse, its mass, outer ra- 
dius, and scale height will take values in the range Md ~ 
1 - IOMq, Rd ~ few hundred AU, and H/R\r^ ~ 0.1 - 0.2, 
respectively. Once the expanding stellar ejecta from the su- 
pernova explosion impacts the disk, a portion of its kinetic 
energy will be thermalized, potentially resulting in bright 
optical or X-ray emission. In this Section, we discuss the 
interaction between the disk and ejecta and its observable 
emission. 
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Figure 4. Prc-explosion properties of the relic star-disk system 
and associated outflows (see Table [TJ. The fractional solid angle 
subtended by the disk is H/R\[i^ ~ 0.1 — 0.2. The inner edge of 
the disk is located at Ra ~ 10 — 30 AU (eq. |25j ) because interior 
to this radius the MRI is fully active and accretion is rapid. The 
outer radius of the disk is Rd ~ 10^ — 10^ AU. The surface density 
between R^^ and Rd obeys E oc R~^ , where the value of 13 1 — 2 
depends on the efficacy of angular momentum transport during 
the irradiation-supported phase ( ^2.3| I. Winds from the star are 
generally fast, with typical speeds ~ lO'^ — 10^ km s~^. Winds 
from the disk, driven either by MHD processes at small radii 
or photo-evaporation coupled with radiation pressure on dust at 
larger radii, are generally slower, with speeds ~ 10 — 100 km s~^. 
Intermediate velocity outflows may originate i n the mixing layer 
betw een the stellar wind and the disk surface l lMatsuvama et al] 
[iooi), or as the result of [shock-mediated] mixing between the 
stellar and disk winds. 



5.1 Ejecta-Disk Interaction 

We begin with a general discussion of the kinematics of 
the disk-ejecta collision. Of the total kinetic energy _Esn ~ 
IqSI _ 10^2 ergs and mass Mcj « M* ~ 10 - IOOMq of the 
SN eject4EI, only a fraction /n = f2d/47r ~ H/R\r^ will in- 
tercept the disk, where fid « A-n{H/R) is the solid angle sub- 
tended by the disk from the star. We consider SN energies 
up to ~ lO^'^ ergs, above the canonical value of 10^^ ergs, 
because VLSNe likely originate from v ery high mass stars 
l|Gal-Yam et al.ll2007l : ISmith et al.ll2007l ') and thei r SNe often 
appear to be intrinsically "hyper-" energetic (e.g. lRest et al.l 
I2OO9). 

In the crude, yet instructive approximation that the in- 
teraction can be treated as a "sticky" collision between two 
masses Mi = fnMi, and M2 = Md, the maximum efficiency 
for converting the SN kinetic energy into thermal energy Eth 
can be shown through conservation of energy and momen- 
tum to be 



-£^th|max 

^max — ~ — J ft 



Mo 



Ml + M2 



H 



rAi + Qd/A 



,(47) 



where Qd > Qo ~ 1 is the Toomre parameter (eq. [T]) eval- 
uated at R = Rd and the constant A — Agi — Ai^ f» 2 — 6.8 



We assume that the SN explosion is largely successful in that 
most of the star is unbound, with only a modest portion of the 
final stellar mass used to form the neutron star or black hole 
compact remnant. 
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Vexp ~ 3000 km/s > W 




Unshocked Stellar / Disk Winds 

(Low Velocity) 



Unshocked Disk 
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Shocked Fast Ejecta 
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Figure 5. Interaction between the outgoing supernova ejecta and 
a relic disk. The bulk of the SN mass and kinetic energy expands 
at a velocity t>exp ~ 3000 km s~^ (eq. 1491 *). The outer, fastest 
portions of the ejecta expand at the speed Dfast ~ 10* km 
and, upon reaching the inner edge of the disk at i? ~ Rim drives 
a strong shock through the midplane. The shock moves through 



at -R ~ -Rin to tjgij ~ 10^ km 



at ij ~ R^, de- 



pending on the disk mass that has been accumulated. If photons 
diffuse out of the midplane and cool the disk faster than the SN 
expansion timescale (viz. tcool ^exp), the light curve directly 
traces the shock-generated thermal energy (eq. 1501 ; Fig. \Ei and 
the SN spectrum may show evidence for CSM interaction, ei- 
ther through intermediate velocity lines (~ 10'^ km s~^) from the 
shocked disk or low velocity lines from the unshocked disk or stel- 
lar/disk winds. If instead tcool 2> tcxp, the energy generated by 
the shocked disk must diffuse through the envelope of the entire 
exploding star (Fig. [Tjl , thereby displaying little direct evidence 
for CSM interaction. 



relates the total disk mass to the surface density at Rd (see 
Appendix A). 

Equation (|47p shows that if the disk does not viscously 
evolve when irradiation-supported (such that Qd ~ Qo ~ 1 
and A — Agi ~ 2 at core collapse; H2.3.ip . then tmax ~ 
H/R\fi^ ~ 10 — 20 per cent. Even if accretion is appreciable 
during the irradiation-supported phase, disks that prefer- 
entially survive have Qd ~ few < A = Ai-^-c — 6.8 at core 
collapse (see Fig. [2|, again resulting in emax ~ 10 per cent. 
If the thermal energy generated by the collision is efficiently 
radiated, the total electromagnetic output may thus reach 
~ Emax-EsN ^ 10^" — 10^^ crgs, significantly higher than that 
from most core-collapse SNe and consistent with the optical 
energies of very luminous SNe. We note that tmax and the 
corresponding luminosity could be even larger if the effective 
solid angle subtended by the disk is larger than our base- 
hne estimate, due to e .g. a puffed up inner disk edge (e.g . 
iDullemond et al.ll200lf) or disk warping (e.g. |Prine:leill996l ). 

We now discuss the radial structure of the disk prior 
to core collapse in order to explore the star-disk interac- 
tion in more detail (see Figure |4] for an illustration). The 
inner edge of the disk occurs at Rin ~ Ra ~ 10 — 30 
AU (eq. because interior to this radius the MRI is 

fully active and accretion is rapid (see i]2.3.2p . Assuming 
that the disk has reached a steady state, the accretion rate 



M oc z/S between Rin and Rd is approximately radially 
constant, where u = ac^/Q oc aTR?^^. If the viscosity is 
very low in the irradiated state and the disk's structure at 
core collapse reflects that of marginal gravitational stability 
( ^2XT1) . then a = agi a T.'^R^^'^ (eq. [6]) and T oc E^i?^ 
(eq. O) and hence E oc i?"^^/^ (see Fig. [XT] in Appendix 
A). On the other hand, if the disk viscously evolves while 
irradiation-supported ( ^2.3.3|) then T oc i?"^/^ (eq. [18]) and 
hence E oc i?^^^''^* for constant a. In general we thus have 
E oc R'l^, where /3 ~ 1 - 2. 

Properly normalized, the pre-supernova surface density 
thus takes the approximate form (for R > Rin) 



E ~ 10 gem" 



Md 



M, 







Rd 



300 AU 



(48) 



where in normalizing we have assumed that Rd ~ lO^iiin 
and neglect the [weak] dependence of the pre-factor on 
/3 ~ 1 — 2. We now scale the disk radius to a somewhat 
larger value Rd = 300 AU, which is more typical of the 
final radius at core collapse (Table [TJ . Exterior to Rd we 
have assumed that E decreases exponentially with radius, 
as would be exp ected from th e isolated viscous evolution of 
an a— disk (e.g. |Pringlelll98ll : see Appendix A). In reality, 
the surface density outside Rd may decrease more rapidly 
with radius if, for instance, the disk is marginally gravita- 
tionally stable ( ^2.3.11 see Fig. IA1[I or if the outer edge has 
been truncated by photo-evaporation (§37TJ or a recent stel- 
lar collision ( >i4.2|l . 

The stellar ejecta expands at an average velocity 



« 3 X lO^kms"^ 



1/2 



Es 



1/2 



1052 ergs J \ lOOM, 







-1/2 



(49) 



and reaches the inner edge of the disk on a timescale 
Rin/vcxp ~ days to a week. Although «exp is the average 
ejecta speed, the expa nsion is homologous with a b road ve- 
locity distribution (e.g. IChevalier Sz Franssonlll994l ). In par- 
ticular, the outer ejecta expands at a typical rate iifaat ~ 10^ 
km s~^ which can be several times larger than Uoxp. Once 
this fast ejecta reaches the inner edge of the disk, a bow 
shock will form and a strong shock will be driven through 
the midplane at the speed Vsh (similar to the interaction with 
giant LBV eruptions; tj4.4|l . This situation is illustrated in 
Figure [S] The shock through the disk generates thermal en- 
ergy at the rate 



{l/2)^Pdvti, X {47V fnR') ~ 7r?/>Ei?< 

, 3 



1.4 X 10''''ergss"V ( 



Vsh 



X 



Md \ / -Rd 
Mq) \ 300 AU 



10* kms-i, 
1^1 exphi?/7?d], (50) 



V Ra 



where pd — E/2// is the midplane density, t/j ~ 1 is an 
efficiency factor, and we have used equation (|48p for E. 

The ratio of the disk density to the mean ejecta density 
Pcj = Mi,/{4iTR^ /3) at radius R can be written as 



pd_ 



1.8Qd\R/Rd 



.(12/7)-l3 



exp[~R/Rd], 



(51) 
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where we have used equation (|19p for H/R. Note that x 
1 is satisfied at all radii for all physical disks (i.e. those 
with Qd ^ 1). Thus, unlike the case of giant LBV eruptions 
described in i|4.4l the efficiency with which the SN ejecta 
shares its energy with the disk material is necessarily high. 

If viscous evolution occurs during the irradiation- 
supported phase (such that P — 15/14 and Qd ~ few) the 
disk mass oc T,R^ oc i?^"^''^'* is concentrated at large radii and 
X remains <^ 1 for i? <C -Rd. The stellar ejecta and the swept- 
up disk mass thus only comparable (and hence deceleration 
becomes significant) once the shock reaches radii ~ Rd- If, on 
the other hand, the disk is marginally gravitationally stable 
at coUapse (such that /3 ~ 27/14 and Qd ~ Qo ~ 1; >i2.3.ip 
appreciable deceleration thus begins almost immediately at 
R ^ -Rin because the disk mass is spread out approximately 
equally per radial decade (see Fig. IA1|) . 



5.2 Two Regimes of Interaction 

In this section we address the emission produced by the disk- 
star interaction. The models we present are not intended to 
describe any supernova in particular and are less detailed 
than other models of CSM interaction available in the lit - 
erature fe.g. IChugai et al.ll2004l : IChugai fc ChevaUeill2006h . 
Rather, our goal is to qualitatively illustrate the diversity of 
observable light curve and spectral behavior that may mani- 
fest from the ejecta-disk interaction and how this may relate 
to observed features of luminous SNe. More detailed radia- 
tion hydrodynamic simulations will ultimately be required 
to confirm and refine the ideas presented here. Although in 
general the ejecta will possess additional sources of energy, 
due to e.g. radioactive decay or residual thermal energy, we 
focus on the luminosity that results solely from the interac- 
tion with the disk. 

Although ijgh (eq. [SD]) represents the power generated 
as the shock passes through the disk, this thermal energy 
may not be immediately radiated. If the disk's opacity is 
K, the vertical optical depth is Td ~ Ek/2 and the photon 
diffusion time out of the midplane (for Rin < R < Rd) is 
given by: 




(52) 

where we have scaled k to the electron scattering opacity 
Kes ~ 0.4 cm^ g~^, a reasonable approximation for the very 
high temperatures T > 10^~* K behind the shock. For E we 
use the post-shock surface density, which is larger than the 
pre-shock value (eq. [48]) by a factor ~ 7, which represents 
the radial compression due to a strong, adiabatic, radiation 
pressure-dominated shocklf] Since the post-shock pressure is 



° Note that if tcool is sufficiently short, the shock will be radiative 
and the post-shock density may become significantly higher than 
just due to adiabatic compression. The adiabatic approximation 
is nevertheless valid on timescales < icooli thereby justifying its 
use in equation II53I I. 
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Figure 6. Light curves from the shock-heated disk in the fast 
cooling case tcool ^ *exp f il5.2.1ll . calculated for a total stellar 
ejecta mass M*, = 30Mq; supernova energy i?sN = 2 X 10^^ 
ergs; and disk mass Md = lAf0 and outer radius Rdj = 300 
AU, respectively. We show calculations for different values of the 
surface density index /9, where S oc (eq. 1481 ). 

dominated by radiation, tdiff.d is also the cooling timescale 
of the disk tcooi (eq. [5]). 

As we show below, the properties of the emission from 
the shocked disk at radius R depends on the ratio of the 
cooling timescale tcooi to the average expansion timescale of 
the bulk of the ejecta tcxp — R/vaxp'- 

^ .02 f^] (_^\ (^^^\ , 
texp ■ \ O.l jyiO^^ergsJ \100Mq J 

(oTifc) (sofxu) (:£) ' ^^^^ 

where we have used equation (|49p for v^xp and have assumed 
5.2.1 Fast Cooling Case 

Low mass and/or radially-extended disks satisfy the condi- 
tion tcool «C iexp- In this case the thermal energy generated 
by the shock escapes the disk before it is engulfed by the ex- 
panding optically-thick star (see Fig. [5]). The shock will thus 
be r adiative and t he swept-up disk will cool into a thin shell 
fe.g. lChugaill200l| j. X-rays generated behind the shock may 
escape t o infinity, producing detectable high-energy e mis- 
sion (e.g. Ilmmler fc Kunt j[2005l : ISchlegel fc Petrdl2006l ). or 
they may be abs orbed and reprocess ed into optical (e.g. 
Ha) emission (e.g. lFransson et al.lll996l ). The observed bolo- 
metric light curve thus directly tracks the shock-generated 
power, i.e. Lboi « E^h. 

Since the swept up disk mass becomes comparable to 
the ejecta mass at radii R < Rd, the cool shell will gener- 
ally decelerate below the expansion speed of the blast wave. 
Photons re-radiated by the disk are thus a possible origin 
for the "intermediate velocity" lines with typical widths of 
a few thousa nd km s~^, as are sometim es observed in Type 
Iln SNe fe.g. lChugai fc Danzigeijll994l '). 

Figure [H] shows the bolometric light curve Lboi = ^^sh in 
the fast-cooling case, calculated from equation (|50p for a disk 
with mass Md = ^Mq and radius Rd = 300 AU; supernova 
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energy Esn — 2 x 10^^ ergs and ejecta mass A/cj ~ 30Mq; 
and for different values of surface density profile index /?. 
The shock heating rate -Esh is calculated using equation (|50p , 
accounting for the deceleration of the s hock speed Vsh(t) us- 
ing th e simple thin shell model given in lChugai fc Danzigerl 
l|l994 ) in which the total momentum Mej /nWexp is conserved. 
We take fn — H/R = 0.1 as the effective scaleheight in this 
one-dimensional model. 

For all values of /3 the total radiated energy is -Bboi ~ 
7 X 10*® ergs or ~ 5 per cent of the SN kinetic energy, 
significantly larger tha n the electromagnetic out put of most 
core-collapse SNe (e.g. iBersten fc Hamuvll2009l ). Also note 
that for P — 0.5, the emission peaks at late times, only 
once the shock reaches the outer edge of the disk at tpcak ~ 
100 days. Such behavior provides one explanation for events 
with long rise tim es until peak emission (e.g. SN 2008iy; 
iMiller et al]|2009l ). 

5.2.2 Slow Cooling Case 

More massive, compact disks instead have tcooi S> texp. In 
this case the shocked thermal energy remains trapped until 
after the disk is engulfed by the expanding, optically-thick 
star. As a result, the shock-generated power iSsh must in- 
stead diffuse out of the entire stellar envelope, which occurs 
on a longer envelope diffusion timescale 



whe re k is the opacity and B ~ 0.07 for a spherical outflow 
(e.g. iPadmanabhan. T.II2OO0I '). 

Figure [7| shows an example model for the bolometric 
light curve in the slow cooling case, calculated for a disk 
with mass Md ~ IOMq and radius Rd — 200 AU; supernova 
energy Ssn = 3 x 10 ergs; ejecta mass Mej = 40Mq . We 
assume the disk density index is /? = 1 and the scaleheight 
H/R = 0.1. We calculate the light curve using a simple one- 
zone diffusion model, in which the radiated luminosity is 
given by 1/ = Eth/tdiS.*, where -Eth is the thermal energy of 
the ejecta (e.g. iLi fc Paczvhskilll99^ ). We evolve Eth in time 
as the star expands, including adiabatic losses due to PdV 
work and time-dependent heating from the shock-generated 
energy Esh. As in the fast-cooling case, we calculate Vahit) 
using a momentum-conserving thin-shell model. We assume 
that the opacity is electron scattering k f» k^s ~ 0.4 cm^ 

The two cases shown in Figure [7] correspond to differ- 
ent assumptions regarding the density profile at the outer 
edge of the disk. In the model shown with a solid line, we 
have assumed that the surface density decreases exponen- 
tially with radius Rd (eq. US]), as expected from the iso- 
lated viscous evolution of an irradiation-supported constant- 
a disk. By contrast, the dotted line shows a model for a 
sharper outer disk edge, as appropriate for a disk which is 
marginally gravitationally-stable ( i]2.3.1l Fig. \Al\i or has its 
outer edge truncated due to disk winds ( i|4.1l i]4.3|) or stellar 
collisions f t|4.2[) . In both cases we find that the total radi- 
ated energy is ~ 3 x 10^° ergs. The radiative efficiency is 
lower than in the fast cooling case because a portion of the 
trapped thermal energy is lost to PdV work before diffus- 
ing out as radiation. Because the ratio t^xp/tdis,* increases 
rapidly with radius (oc R^), adiabatic losses are, however, 



Slow Cooling Case (t,„„| » t.,„) 
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Figure 7. Light curve powered by the shocked disk in the slow 
coohng case tcool ^ 'cxp ( ^5.2.2[ . calculated for a total stellar 
ejecta mass M* = 40A/q; supernova energy -BsN = 3 x 10^^ ergs; 
and disk mass Md = IOMq and outer radius Rd = 100 AU, re- 
spectively. We take /3 = 1 in the disk density profile. The model 
shown with a solid line assumes that the disk surface density out- 
side Rd decreases exponentially with radius, as expected from the 
isolated viscous evolution of an irradiation-supported, constant-a 
disk (eq. 1481 ). The calculation shown with a dashed line instead 
assumes a sharper outer disk edge, as would be expected if the 
disk resides in a state of marginal gravitational instability at col- 
lapse ( il2.3.1l see Fig. lAlj l or if the outer edge has been truncated 
due to, e.g., photo-evaporation ([|4!TJ or a recent stellar collision 

(EM- 

mild compared to those incurred from the compact stellar 
surface. 

Although we have discussed the fast and slow cooling 
cases as if they are mutually exclusive, in general tcooi/ioxp 
will vary with radius and the resulting emission will be more 
complex than the two limiting cases described above. For 
instance, tcooi/tcxp > 1 may be satisfied at small radii, in 
which case a portion of the thermal energy generated by the 
shocked inner disk will be trapped behind the stellar photo- 
sphere. However, if fcooi/toxp oc R~'^ ;^ 1 is satisfied once as 
the shock reaches the outer disk, a significant portion of the 
thermal energy could also be radiated promptly. This type 
of beh avior could e x plain "transition" events like SN 2006tf, 
which ISmith et all pOOS') argue evolves from an optically- 
thick, "diffusive" phase to emission which shows more direct 
evidence for CSM interaction at later times. 



6 DISCUSSION 

6.1 Long-Lived Disks around Very Massive Stars 

A straightforward prediction of the relic disk model is 
that a small fraction of very massive stars should possess 
a proto-stellar disk at radii of tens to hundreds of AU 
throughout their entire lifetime (including putative post- 
main sequence LBV and Wolf-Rayet phases). The chief 
difficulties in testing this prediction are that (1) massive 
stars are typically far away, and the ir environments tend 
to be messy, crowded, a nd obscured (jCesaroni et al ] |jiHd : 
IZinnecker fc Yorkd |2007| ): and (2) because very luminous 
SNe are very rare, only a small fraction of massive stars 
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must retain a disk in order to account for VLSNe via 
CSM interaction with a relic disk. For instance, if a frac- 
tion / ~ 10"^ - 10"* of all core collapse events pro- 
duce VLSNe, then at any time we would expect to find 
N ~ /i?sNtiife ~ 3 — 300 progenitors in a galaxy like the 
Milky Way, where -Rsn ~ 10~^ yr~^ is the total core col- 
lapse rate and tutc ~ 3 Myr. We now discuss whether such 
progenitor systems actually exist. 

It is clear that long-lived disks are not present 
around the majority of [what are normally considered] 
main sequence O stars. The emission line profiles of some 
main sequence O stars (the "Oe" and "Onfp" classes) 
have long been interpreted as indicating the presence of 
an ionized disk [C onti & Lccp 1974). Spex^trogolarimetry 
of O stars (iHarries et all l2002l: [Vink et all 120091 '). LBVs 



(e.g. Schulte-Ladbeck et al.l [1993), and Wolf-Rayet stars 



(e.g. IVin'n 2007 ') also suggest that the environments of mas- 



sive stars are asymmetric in some cases. These probes are, 
however, primarily sensitive to scales of only a few stellar 
radii, and the observed polarization may be caused by intrin- 
sic asy mmetries in the stellar wind, p erhaps due to rapid ro- 
tation (|Biorkman fc Cassinelliiri993h . Furthermore, massive 
irradiated disks around otherwise unobscured stars would 
likely produce extremely bright emission line (e.g. Ha) lu- 
minosities, which are not observed. 

A long-lived disk should also be a strong source 
of infrared (IR) emission. Since the disk scaleheight is 
H/R ~ 0.1, if the disk is flared then at least ~ 10 per 
cent of the stellar luminosity will be absorbed and re- 
radiated as thermal blackbody and d ust emission peaked 
at near- to mid-IR wavelength s (e.g. IChiang fc GoldreichI 
Il997l : iRafikov fc De Collel mm . Excess IR emission is in 
fact common from very massive stars during es s ential ly all 
stages of evolution. For instance, iBarniske et al.l (|2008l ) find 
evidence for mid-IR dust emission down to scales < 100 AU 
around two nitrogen-rich WN stars near the Galactic cen- 
ter, which they infer to be a mong the most lumino us stars 
in the Milky Way (see also iRaiagopal et all 120091 '). Wolf- 
Rayet (WR) stars, especially of the carbon sequence (WC), 
commonly radiate a significant fracti on (~ 0.1 — 10%) of 
their total luminosit y in the IR (e.g. IWilliams et al.|[l987l : 
iHadfield et al] 120071 '). In cases when the infrared emission 
is p artially resolved, the inferred radii are ~ 30 — 300 
AU (|Dvck et al.l ll984:'Ragl and fc Rich ichi''l999'). consistent 
with the expected scale of a relic disk. Excess IR emission is 
also common from LBVs and, when margin ally resolved, size 
scale s of tens of AU are commonly inferred l|Raiagopal et al] 
l2007l ). Although these hints are suggestive, the mid-IR emis- 
sion from massive stars is generally well-modeled as originat- 
ing from a quasi-sph erical, dusty outflow, without the nee d 
to invoke a disk (e.g. IWilliams et al.|[T987l : ICrowtheill2003l '). 

Although there appears to be no evidence for long- 
lived disks around unobscured stars, it is important to keep 
in mind that the time-sequence and durations of the ob- 
serva ble phases of massive stellar evolution remain uncer- 
tain (|Zinnecker fc Yorkc 200?! ). When massive stars are very 
young and still actively accreting, they may be observed as 
"hot molecular cores," with large mas ses of warm and dense 
gas (|Kurtz et al.ll2000l : [c esaronil l2005l l. Once external accre- 
tion shuts off, howeve r, the system is thought to transition to 
a compact HII phase (|Kurt j2005l : lHoare et al.ll2007i '). during 



which the central star begins ionizing pockets of surrounding 
gas. 

After the HII phase, the system is generally thought 
to emerge as an unobscured main sequence star. For mas- 
sive stars, however, the HII phase is known to last a sig- 
niflcant fract ion (~ 10 per cent, on av erage) of the stel- 
lar lifetime (Wood & Churchwcll 1989). In order to ex- 
plain the anomalously lo ng lifetimes of compact HII regions, 
iHollenbach et all (|l994l ) proposed that they are continually 
replenished with material from a photo-evaporating disk; 
this model in fact re mains a leading explanation for hyper- 
compact HII regions (|Ketcll2007l : iNielbock et al.ll2007l ') with 
electron number densities rie > 10^~® cm~^ and inferred 
sizes < 10^ AU (e.g. lKurtz et al.ll2000l : [be Free et al.ll2005l ). 
Although disks around the most massive stars may not 
completely photo-evaporate (§33} or be stripped by stel- 
lar winds ( i]4.3|) . the mass loss rate from the disk can still 
be quite high, such that several solar masses could be lost 
from the disk throughout the stellar lifetime. By producing a 
dusty shell around the star, a continuous disk outflow could 
in principle cause even a ~ Myr old star to appear to be in 
a younger stage of evolution. 

To summarize, if VLSNe indeed originate from relic 
proto-stellar disks then their progenitors should manifest as 
conspicuous sources of IR emission and possibly as long-lived 
hyper-compact HII regions. If unobscured, such systems 
should be exceedingly bright sources of optical (e.g. Ha) 
line emission. Although no such progenitor population has 
yet been identified, only a small fraction of massive stars 
produce VLSNe and evolved stars with relic disk systems 
could in principle be confused with younger embedded sys- 
tems. We note that only quit e recently are blind samples of 
hyper-compact HII regions (|Murphv et al.l '2010'l. a nd IR- 
selected samples of evolved massive star s (Hadfie ld et al.l 
l2007l : IShara et all l2009l : ICarev et al.l |2009| ) and the general 
popu l ation of self-obscu red massive stars l|Thompson et al.l 
I2OO9I : iKhan et al.ll2010l ). being assembled. The possibility 
that a population of relic disk systems may exist without 
our knowledge is underscored by the recent discovery of a 
"hidden" population of massive stars with the Spitzer Space 
Telescope (jWachter et al.ll2010l '). 



6.2 Very Luminous Supernovae 



Many VLSNe a re classified as Type Iln (|Schlegell Il990l : 
lFilippenkolll997l ). which are characterized by (1) slowly- 
deca ying light curve s that can persist for decades 
(e.g. lMiller etaLll2010l ): (2) very narrow (widths of tens to a 
few hundred km s~^) H emission lines, which are thought to 
arise from photo-ionized unshocked CSM; (3) intermediate- 
width (typically 1 — 5 x 10^ km s~ ^) H lines, thought to aris e 
from shock-accelerated CSM (e.g. lChugai fc Danzigeijll994l ). 
(4) and, sometimes, luminou s X-ray and radio emis sion for 
years after the explosion (e.g. Ilmmler fc Kunt jlioosi ). These 
characteristics are in addition to the broad lines with widths 
~ 10* km s~^ which are observed in most SNe. The progen- 
itors oOIn_SNe_axe_£rotoblj;^ the most massive stars 
(e.g. iGal-Yam fc Leonard! [2OO9I ). 

As described in iJH many of the properties of the most 
luminous SN Iln can be produced via the interaction of the 
supernova shock with a massive, relic proto-stellar disk (see 
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Figure [5]). Broad lines may originate from the shocked fast 
supernova ejecta, which occupies the majority (~ 80 — 90%) 
of the solid angle around the exploding star; they may be 
powered by either the fast ejecta interacting directly with 
the dense wind from the star or disk wind; or through the 
absorption of X-rays radiated at mid latitudes as the shock 
passes through the disk. Intermediate lines may originate 
from the shocked disk, which is initially accelerated by the 
fast ejecta, but decelerates once it has swept up apprecia- 
ble mass. Since the mass of the disk and the sector of the 
star which intercepts the disk are approximately equal (see 
eq. [47)). intermediate lines with widths 1/2 of the bulk 
ejecta speed Uexp (eq. [49]) around the time of peak emission 
are thus a natural prediction of the relic disk model. 

The narrow line emission in SNe Tin is generally well- 
modeled as originating from the outer, unshocked portion 
of the same outflow which is responsib le for the bulk of 
the luminosity (e.g. IChugai et al.l [20o3 ). However, in the 
relic disk model, narrow line emission must originate from 
a region distinct from the disk itself. Disk outflow due to 
photo-evaporation ( H4.ip or stellar wind stripping ( H4.3P are 
generally too slow to explain the line widths of up to sev- 
eral hundred km s~^ that are observed in luminous SN Iln. 
However, in principle the faster stellar wind and the slower 
disk material could mix, resulting in an outflow with in- 
termediate velocities. Indeed, empirically broad radio re- 
combination lines with widths 30 — 200 km s~^ are ob- 
served in hyper-cornpact HII regions dCaume et all 1 19951 : 



jrn pa 

IShepherd et'Hlll995l : ls"ewilo et al.ll2004l ). which are thought 



to result from p hoto-evaporating disks (|Hollenbach et ahl 
ll994l : lKetdl2007f ). 

The majority of SN Iln probably cannot result from 
relic disks. Since SN Iln represent up to ~ 10 per cent of core 
collapse SNe, the implied population of progenitors would 
probably have been identified ( i]6.ip . On the other hand, SN 
Iln are an inhomogeneous population and the origin of the 
GSM in the general population may be different than in the 
most luminous cases. 

We also emphasize that not all VLSNe are Type 
Iln. Many of t he mo st energetic eve nts such as 2 005ap 



2007f) and 2008es l|Miller et al.l l2009l : 

iGezari et al.l 20091 ) lack intermediate width lines (~ 2000 



( Quimbv et al 



km s~^), while the high-re dshift PTF events show no evi- 
dence for H or He at all (Quimby et al.i! 2009i ). Although on 
the surface these observations appear difficult to explain via 
interaction with a hydrogen-rich proto-stellar disk, evidence 
for GSM interaction could be masked if the disk is particu- 
larly massive and/or compact because it may not have time 
to radiate intermediate-width lines before being engulfed by 
the expanding ejecta (see il5.2.2l and Fig.[7l). Indeed, our cal- 
culations in 33] suggest that massive, compact disks which 
satisfy tcooi S> tcxp (eq. [53]) are only expected to survive 
around the highest mass (and, hence, rarest) stars which 
may have the most energetic explosions. 

6.3 "Hybrid" Type I/IIn SNe 

SN 2002ic (|Hamuv et all 120031 ') was classified as a pecu- 
liar Type la SN, but with the unusual feature of narrow 
and broad components of Ha emission, suggesting the pres- 
ence of strong H-rich GSM interaction as in SNe Iln. This 
"hybrid" Type la/IIn SN was interpreted bv iHamuv et al.l 



(I2OO3I ) as resulting from the thermonuclear disruption of a 
white dwarf into a hydrogen-rich envir onment, perhaps re - 
sulting from a pre-explosion A GB wind (iHarnuv et al.ll2003l ) 
or common envelope ejection (jLivio fc RiessI 20031 ). The re- 
quired GSM mass to explain the observed light curve was 
up to several solar masse s on a radial scale ~ 500 AU 
(iGhugai fc Yun gclson'2004'). SN 2005gj (jPrieto et al.ll2005l : 
lAldering et al. |2 006 ) also showed hybrid Type la/IIn char- 
acteristics similar to 2002ic, which again appeared to require 
a significant quantity of mass on a radial scale of a few hun- 
dred AIL 

'B enetti et al ] Jlool) recently questioned the identifica- 
tion of 2002ic as an intrinsic Type la SN, suggesting it should 
instead be classified Type Ic, events which are thought to 
result from the core collapse of WR stars. Although it is 
unclear why a star would lose its entire hydrogen enve- 
lope just years prior to core collapse, hydrogen-rich ma- 
terial around an evolved star at a few hundred AU could 
be explained as a relic proto-stellar disk. Spectropolarime- 
try of 2002ic supports this hypothesis by showing that the 
geometry was likely aspherical and, possibly, equ atorially- 
concentrated (iKotak et alj|2004l: IWang et aLllioM). Analv- 
sis of the light curve ( Wood-Vasev et al.ll2004l : Uenishi et al.l 
l2004f ) further showed a delay between the time of explosion 
and the onset of the GSM interaction, suggesting the pres- 
ence of a "cavity" between the star and the GSM. A lower- 
density region in proto-stellar disks is expected interior to 
the MRI-active radius Ra ~ 30 AU (eq. [55]; see Fig. [4]). 



7 CONCLUSIONS 

Gircumstellar interaction is a promising mechanism for pow- 
ering some of the most luminous and energetic SNe yet dis- 
covered. However, if the requisite GSM is the result of a 
pre-SN stellar eruption, up to several solar masses of mate- 
rial must be ejected just years— decades prior to explosion. 
Although in principle giant LBV-like eruptions or pair insta- 
bilities are capable of producing such prodigious mass loss, 
there is currently no explanation for why the eruption and 
the stellar explosion are so nearly synchronized. 

In this paper we have developed an alternative hypothe- 
sis that the bulk of the GSM is not an outflow (thus avoiding 
the "coincidence" problem), but rather a pre-existing, long- 
lived disk left over from stellar birth. One virtue of our model 
is that a GSM mass up to ~ IOMq is explained as the max- 
imum disk mass that can be supported against self-gravity 
by stellar irradiation (eq. [20]). The GSM radu ~ 10^ - 10^ 
AU required by observations are also explained naturally 
by (1) the maximum outer disk radius (~ 100 — 200 AU) 
that is stable to fragment ation during the embedded phase 
(|Kratter fc Matzner|[200^ : gUTJ; and (2) the limited subse- 
quent viscous evolution experienced by disks around massive 
stars ( ^2.31 !j3]). Indeed, we find that massive disks may pref- 
erentially survive until core collapse around the most mas- 
sive stars ( i]4.5[) for the following combination of reasons (see 
also Fig. [3]): 



• Massive stars live short lives, thereby allowing less time 
for accretion and other disk dispersal processes to act. 

• Massive stars begin their irradiation-supported phase of 
evolution with more massive disks (eq. [20]), which generally 
take longer to disperse. 
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• Disks around massive stars have larger surface densities 
at the onset of irradiation-supported accretion (eq. [22]). 
This shields the midplane from external sources of ion- 
ization, resulting in extensive and sustained dead-zones in 
which the magneto-rotational instability is quenched and 
angular momentum transport is inefficient. The powerful 
stellar winds from massive stars appear likely to exclude 
interstellar cosmic rays, leaving X-rays (which penetrate to 
a shallower depth) as the chief source of surface ionization 
(eq. [23]). 

• Due to their slower expected viscous spreading and 
deeper gravitational potential, disks around massive stars 
are less likely to be dispersed via photo-evaporation ( i]4.ip . 

• Massive disks can be denser than the ejecta from giant 
LBV eruptions, suggesting they are more likely to remain 
bound to the star ( ^4.4p . 

Our conclusion that the probability of disk survival rapidly 
increases with stellar mass is consistent with evidence that 
very luminous Iln SNe originate from the most massive stars 
l|Gal-Yam et al.ll2007l : fSmith et al.ll2007l '). We do not, how- 
ever, suggest that the disk survival fraction approaches unity 
in any stellar mass range, as external influences, such as the 
presence of a massive binary companion or a particularly 
dense stellar environment, may preclude disk survival in the 
majority of cases. 

Several of the observed properties of VLSNe can be ex- 
plained by the collision between the supernova ejecta and a 
relic disk (see Figs.|4]and[5]). The nominal radiative efficiency 
of the disk-ejecta interaction is ~ 10 per cent (eq. [13), re- 
sulting in a total electromagnetic output ~ 10^" — 10^^ ergs 
for a SN with kinetic energy ~ 10^^ — 10^^ ergs, substan- 
tially more optically-energetic than a normal Type II SN. 
For lower mass or radially-extended disks, the light curve 
directly traces the shock-generated thermal energy ( i]5.2.1l 
Fig. [ni) plausibly resulting in copious "intermediate" width 
line emission, as observed in SN Iln. For more massive, com- 
pact disks the shock-generated thermal energy must instead 
diffuse out of the entire stellar envelope ( tj5.2.2l Fig. [7|), 
which may explain the most luminous events (e.g. 2005ap), 
which show little or no direct evidence for GSM interaction. 

Although there is no observational evidence for long- 
lived disks around massive stars, distinguishing between 
a truly young stellar object proto-star and an older sys- 
tem with a relic disk may be nontrivial. Ultimately test- 
ing whether bona fide disks exist around evolved stars in 
sufficient numbers to explain VLSNe will require acquir- 
ing unbiased samples of embedded evolved stars and hyper- 
compact HII regions, distinguishing evolved from youthful 
systems, and, eventually, pinning down definitive disk or 
outflow signatures with e.g. atomic {Ercolano & Owen 2010) 
and molecular line diagnostics (e.g. Kr umholz et al.l 120071 : 
ICarr fc Naiit^l2008f ). Though challenging, such a feat will 
be aided by the improved capabilities of present and upcom- 
ing observatories such as the EVLA, ALMA, Herschel, and 
the James Webb Space Telescope. 
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APPENDIX A: CALIBRATION OF THE RING 
MODEL 

This Appendix describes how the constants A, B, and / 
employed in the ring model in are determined. Most of 
this material is from Metzger, Piro, & Quataert (2008), but 
we repeat it here for convenience. 

The surface density E of an axisymmetric disk in a Ke- 
plerian potential with constant total angular momentum 
evolves according to a diffusion equation (e.g., Frank et 
al. 2002): 



9E 
dt 



3d_ 

r dr 



1/2 d I i/2\ 



(Al) 



where v is the kinematic viscosity. Assuming that v depends 
only on radius as a power law, viz: v = foir / Rd,o)" , equation 
(|A1|) is linear and, for an initial surface density distribution 
E(r, f = 0) = {AId,o /2-KRdfl)S{r — Rd,o) which is narrowly 
peaked about the radius Rdfi, the solution (for n < 2) is 
given by 



S(r,t) = 

Md,o(l-n/2) 
ttRI (,j;("+i/4)r 



exp 





'2x^-^1^' 






T 



,(A2) 



where Mdfi is the initial disk mass, x = r/Rdfi, r = 
f[12!-'o(l — n/2Y / R^ and is a modified Bessel func- 
tion of order m. For small argument y ^ 1, Im{y) takes 
the asymptotic form Im — {y /2)"^ /r{m + 1), where F is the 
Gamma function; thus, for late times or small radii such that 



r > 2x 



l-n/2 



, equation (IA2p reduces to 



Md,o {l-n/2) 



Ff^ 



2n 1 
2nJ 



■ exp 



-{1 + x^-^) 



(A3) 



Most of the mass in the disk is located near the radius where 
the local mass Md oc Er^ peaks; using equation (|A3|I , at late 
times this radius is found to be rpoak = -Rd.o''"^'^'^ . Hence, 
equation (|A3|l becomes valid near rpeak for r S> 1. 

The constant A, which relates the total disk mass at late 
times from the exact solution of equation (|A1|) to the mass 
defined by 7rE(rpcak)'"pcaki can be calculated from equation 
(lA3ll to be 



A{t > 1) 



/■ cx 




2n'Erdr 



7rS(rpeak)r 



pealc 



2e 



(A4) 



Similarly, the constant B, which relates the total disk angu- 
lar momentum at late times from the exact solution to that 
estimated by ■KT,rp^^^^{GMrpca.k)^^'^ , is given by 
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B(r > 1) 



Jo 



7rE(rpcak)rp{f^ 



2e 



-5 - 2n- 
A - 2n. 



(A5) 



From mass continuity, the radial velocity is given by 
-3 d f ^ 1/2] _ — Si'o 1 d 



T,r^/2 dr 



[Ex"+^/2] , (A6) 



which, using equation 
radii 



Rd,o Ea;i/2 dx 
3]) , gives the accretion rate at small 



Min = -27rEriIr|^^2;cl-"/2 

Md,o 3(1 - n/2) 

^1.0/^0 r[(5 - 2n)/(4 - 2n)] 



exp[— l/r]r 



(A7) 

Equation (|A7|I is easily checked by noting that Mindt = 
Md,o, which shows that the entire initial disk eventually ac- 
cretes onto the central object. In tj3.1l we introduced the 
following prescription for evolving the disk mass: 



Mi = 



(A8) 



where, in terms of the viscosity prescription adopted above, 
tvisc = R^d/v = tvisc,o(-Rd/-Rd,o)^~" and tvisc.o = Rdfil^o is 

the initial viscous time. Assuming that the total disk an- 

1/2 

gular momentum remains constant, viz. J cx: MdfiRdQ = 
{B / A)MdB^J^ the solution to equation HA8|) is given by 

Md{t) = Md,o[l+(4-2n)(BM)(^-2"'/(iAvisc.o)]-'''''"'"'(A9) 

In our evolutionary calculations we set / so that the ac- 
cretion rate from the exact solution to equation (|A1[) (A/in; 
eq. I A 71 1 matches the solution to equation ([ASP at late times 
(i.e., in the self-similar limit). This requires 

/ = 3(1 - n/2) (AlO) 



Al Irradiation-Supported Accretion 

For an irradiation-supported, constant-a disk v oc H^Q, oc 
r^^'^^ (eq. [19]); thus, n = 15/14, / = /;„ ^ 1.4, A = Ai„ ^ 



6.80, and B = Bir 



5.85. 



A2 Gravito-turbulent Accretion 

For a gravito-turbulent disk u = a^^CsH oc E^r^^ 
Since v in this case depends explicitly on E as well as r, 
equation ([Aip is a no longer linear and the analytic meth- 
ods employed above do not strictly apply. In Figure we 
show E(r, t), calculated by solving equation (|A1|) numeri- 
cally for a disk with viscosity v cx E^r^^ and a mass Mdfi 
initially concentrated at the radius Rdfi- We show E(r) at 
several times, normalized to the initial viscous time at -Rd,o- 
Figure lAll shows that as the disk spreads from its initial 
configuration, a density profile E oc r"^^''^, corresponding 
to a radially-constant accretion rate M oc i/E, is rapidly 
established between the inner and outer edge of the disk. 
Note that the disk mass (oc Er'^ oc r~^^^) is concentrated 
approximately equally per unit decade in radius. This im- 
plies that the constant A (eq. [A4j ') is well approximated as 
^gi = Md/-KRl « ln(i?d/i?in), where i?in ^ Ra - 2>Q 

AU is the inner edge of the disk (eq. [25]). For a typical 




r/Rd,o 

Figure Al. Time evolution of a gravito-turbulent disk with vis- 
cosity V oc Ti^r^'', calculated by solving equation I)A1|| numeri- 
cally. The disk is is initially concentrated at the radius r = Rd , 
and snapshots of the radial profile of the local disk mass Trr^S 
are shown at various times in units of the initial viscous time. 
Note that rapidly into the disk's viscous evolution, the disk edges 
become sharp and mass is spread out approximately equally per 
decade in radius. 



value Rd ~ lO-Rin, we thus estimate that A = Agi ~ 2 for 
gravito-turbulent disks. 

Although the disk evolution is not strictly described 
by the analytic formalism applied above to the irradiation- 
supported case, we may still estimate the value of the con- 
stant /gi by recognizing that near the peak radius Rd the sur- 
face density locally obeys E oc Md/Rd, implying that (near 
the peak) we have ly oc E^r^"" oc M^r"^ oc for constant 



Jd oc MdR]J^. Forn = equation (fSlO)) gives / = /gi = 3, 
the value which we thus employ in our evolutionary calcu- 
lations in ij3] 

We note, however, that it is not clear a prion that the 
gravito-turbulent disk evolution is indeed dominated by the 
viscosity at radii ~ Rd- The edges of the disk profile in Fig- 
ure lAll are very sharp because the viscous spreading time 
oc /v cx E~^ becomes large in regions where E becomes 
low. Low viscosity near the disk edge also causes the disk 
to spread slower than the prediction Rd oc t^^'^ of the self- 
similar model (see Appendix B), which assumes that the 
disk evolution is dominated by regions of the disk where the 
majority of mass and angular momentum reside. This slug- 
gish evolution is, however, probably unphysical. For one, the 
assumption of a gravito-turbulent viscosity becomes invalid 
near the disk edge in regions where the temperature de- 
creases below the floor set by external radiation (see eq. [20] 
and surrounding discussion). Furthermore, even the assump- 
tions inherent in a height-integrated accretion model become 
questionable regions where the surface density decreases so 
sharply that radial gradients exceed vertical gradients. Thus, 
although the qualitative features in Figure [S] such as the fiat 
mass profile and sharp disk edges are plausibly physical, 
the sluggish evolution is likely an artifact of idealized as- 
sumptions. We conclude that a self-similar model using tvisc 
evaluate d at _R Rd m ay thus remain a reasonable approx- 
imation (|Pringlelll99ll '). 
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APPENDIX B: ANALYTIC SELF-SIMILAR 
SOLUTIONS 

Neglecting details such as stellar mass-loss, the late-time 
evolution of the disk properties from our calculations in ^ 
asymptote to analytic self-similar power-law solutions. 

From equation (|A9|) we observe that at times t ^ fvisc.o 
the disk mass and radius evolve as 

-l/(4-2n) 



Md/Md,o = Cm 



t 



^visc.O 



where 
Cm ~ r 
and 

Rd/Rdfi 

where 

Cr = : 



5 - 2n 



4 - 2n 



((3/4)(4-2n) 



2\-l/(4-2n) 



Md A 
Md,o B 



= Cr 



t 



1/(2-") 



5 - 2n 



2n 



((3/4)(4-2n) 



2\l/(2-n) 



(Bl) 



(B2) 



(B3) 



(B4) 



and equation (IB3|I follows from angular momentum conser- 
vation. 

For an irradiation-supported disk n = 15/14, Md oc 
t'"^^^^ and iid CK t^"*/^^. Although gravito-turbulent disks 
obey a non-linear diffusion equation, we find n = Q us- 
ing disk properties evaluated locally at radii ^ Rd (see 
Appendix A). In this case we thus have Md cx and 
-Rd oc t^^'^ . As described in ^ all of the other disk properties 
a± R = Rd follow from Md{t) and Rd{t). The full set of self- 
similar solutions for the gravito-turbulent and irradiation- 
supported a-disk are given in equations and [28|-|34). 
respectively. 
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